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Abstract 

We present an overview of recent experimental and theoretical ad- 
vances in our understanding of the spin structure of protons and neutrons. 
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1 Introduction 

Attempting to understand the origin of the intrinsic spin of the proton and 
neutron has been an active area of both experimental and theoretical research for 
the past twenty years. With the confirmation that the proton and neutron were 
not elementary particles, physicists were challenged with the task of explaining 
the nucleon's spin in terms of its constituents. In a simple constituent picture 
one can decompose the nucleon's spin as 

Jf = + + + Lf = i . (1) 

where Sz and Lz represent the intrinsic and orbital angular momentum respec- 
tively for quarks and gluons. A simple non-relativistic quark model (as described 
below) gives directly 5"! = ^ and all the other components = 0. 



2 



Because the structure of the nucleon is governed by the strong interaction, 
the components of the nucleon's spin must in principle be calculable from the 
fundamental theory: Quantum Chromodynamics (QCD). However since the 
spin is a low energy property, direct calculations with non-perturbative QCD 
are only possible at present with primitive lattice simulations. The fact that 
the nucleon spin composition can be measured directly from experiments has 
created an important frontier in hadron physics phenomenology and has had 
crucial impact on our basic knowledge of the internal structure of the nucleon. 

This paper summarizes the status of our experimental and theoretical under- 
standing of the nucleon's spin structure. We begin with a simplified discussion 
of nucleon spin structure and how it can be accessed through polarized Deep 
Inelastic Scattering (DIS). This is followed by a theoretical overview of spin 
structure in terms of QCD. The experimental program is then reviewed where 
we discuss the vastly different techniques being applied in order to limit possi- 
ble systematic errors in the measurements. We then address the variety of spin 
distributions associated with the nucleon: the total quark helicity distribution 
AS extracted from inclusive scattering, the individual quark helicity distribu- 
tions (flavor separation) determined by semi-inclusive scattering, and the gluon 
helicity distribution accessed by a variety of probes. We also discuss some ad- 
ditional distributions that have have recently been discussed theoretically but 
are only just being accessed experimentally: the transversity distribution and 
the off-forward distributions. Lastly we review a few topics closely related to 
the spin structure of the nucleon. 

A number of reviews of nucleon spin structure have been published. Fol- 



lowing the pioneering review of the field by Hughes and Kuti |177| which set 



the stage for the very rapid development over the last fifteen years, a number 



of reviews have summarized the recent developments ||7^, 216, 176, 105|. 
Also Ref. presents a detailed review of the potential contribution of the 
Relativistic Heavy Ion Collider (RHIC) to field of nucleon spin structure. 

1.1 A Simple Model for Proton Spin 

A simple non-relativistic wave function for the proton comprising only the va- 
lence up and down quarks can be written as 

\p T) = -^{2\u UUl)-\uUld^)-\uiuU';)- (2) 

where we have suppressed the color indices and permutations for simplicity 
but enforced the normalization. Here the up and down quarks give all of the 
proton's spin. The contribution of the u and d quarks to the proton's spin 
can be determined by the use of the following matrix element and projection 
operator: 

= {pUOu^\p^) (3) 

aT = ^(i + ^3)(i + ^3) . (4) 
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where the matrix element gives the number of up quarks polarized along the 
direction of the proton's polarization. With the above matrix element and a 
similar one for the down quarks, the quark spin contributions can be defined as 

Au^v) -u^ = i((a3) + (asTs)) - ^ , (5) 
M = d^ -d^ = \{{a:,)-{a,T,)) = -\. (6) 

Thus the fraction of the proton's spin carried by quarks in this simple model is 

AE = Au + Ad + As = ((73) = 2 Jf = 1 , (7) 

and all of the spin is carried by the quarks. Note however that this simple 
model overestimates another property of the nucleon, namely the axial-vector 
weak coupling constant qa- In fact this model gives 

^ = (a3r3)=Au-Ad=^ , (8) 
gv 3 

compared to the experimentally measured value of QaI gv = 1.267± 0.004. The 
difference between the simple non-relativistic model and the data is often at- 
tributed to relativistic effects. This "quenching" factor of ~ 0.75 can be applied 
to the spin carried by quarks to give the following "relativistic" constituent 
quark model predictions: 

AS « 0.75 , 
Au« 1.0, 
Ad« -0.25, 

As « . (9) 



1.2 Lepton Scattering as a Probe of Spin Structure 

Deep-inelastic scattering (DIS) with charged lepton beams has been the key 
tool for probing the structure of the nucleon. With "polarized beams and targets 
the spin structure of the nucleon becomes accessible. Information from neutral 
lepton scattering (neutrinos) is complementary to that from charged leptons but 
is generally of lower statistical quality. 

The access to nucleon structure through lepton scattering can best be seen 
within the Quark-Parton Model (QPM). An example of a deep-inelastic scat- 
tering process is shown in Fig. |^. In this picture a virtual photon of four- 
momentum g^(with energy v and four- momentum transfer = —q^) strikes 
an asymptotically free quark in the nucleon. We are interested in the deep- 
inelastic (Bjorken) limit in which and v are large, while the Bjorken scaling 
variable xb — /2Mv is kept fixed (M is the nucleon mass). For unpolarized 
scattering the quark "momentum" distributions — qi{x) = u{x), d{x), s{x), ... — 
are probed in this reaction, where x = xb is the quark's momentum fraction. 
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Figure 1: Deep Inelastic Scattering in the Quark-Parton Model 

From the cross section for this process, the structure function Fi{x) can be 
extracted. In the quark-parton model this structure function is related to the 
unpolarized quark distributions via 

i 

where the sum is over both quark and anti-quark flavors. With polarized beams 
and targets the quark spin distributions can be probed. This sensitivity re- 
sults from the requirement that the quark's spin be anti-parallel to the virtual 
photon's spin in order for the quark to absorb the virtual photon. With the 
assumption of nearly massless and collinear quarks, angular momentum would 
not be conserved if the quark absorbs a photon when its spin is parallel to the 
photon's spin. Thus measurements of the spin-dependent cross section allow the 
extraction of the spin-dependent structure function gi (x) . Again in the quark- 
parton model this structure function is related to the quark spin distributions 
via 

i 

The structure function gi is extracted from the measured asymmetries of the 
scattering cross section as the beam or target spin is reversed. These asymme- 
tries are measured with longitudinally polarized beams and longitudinally (^||) 
and transversely {A±) polarized targets (see Sec. ||). 

Beyond the QPM, QCD introduces a momentum scale (Q^) dependence 
into the structure functions (eg. Fi{x,Q'^) and gi{x,Q'^)). The calculation of 
this dependence is based on the Operator Product Expansion (OPE) and 
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the renormalization group equations (see eg. Ref. [261, 18C, ^). We will not 
discuss this in detail, but we will use some elements of the expansion. In par- 
ticular, the expansion can be written in terms of "twist" which is the difference 
between the dimension and the spin of the operators that form the basis for 
the expansion. The matrix elements of these operators cannot be calculated in 
perturbative QCD, but the corresponding Q^-dependent coefficients are calcu- 
lable. The lowest order coefhcients (twist-two) remain finite as ^ oo while 
the higher- twist coefficients vanish as ^ oo (due to their 1/Q^ dependence). 
Therefore, the full dependence includes both QCD radiative corrections (cal- 
culated to next-to-leading-order (NLO) at pres ent) and higher-twist corrections. 



The NLO corrections will be discussed in Sec. 3.5 



1.3 Theoretical Introduction 

To go beyond the simple picture of nucleon spin structure discussed above we 
must address the spin structure within the context of QCD. We discuss several 
of these issues in the following sections. 



1.3.1 Quark Helicity Distributions and gi{x,Q'^) 

In polarized DIS, the antisymmetric part of the nucleon tensor is measured, 

^^"^hj e^^'A(^5|[J^(0, J.(0)|P5) , (12) 

where \PS) is the ground state of the nucleon with four-momentum P^^ and 
polarization {P ■ S = 0), and J'^ is the electromagnetic current. The anti- 
symmetric part can be expressed in terms of two invariant structure functions, 

^ _^eM-/3^^ [Gi(z., Q^)Sp/M^ + G2{<^, Q^){SpvM - PpS ■ q)/M*] . 

(13) 

In the Bjorken limit, we obtain two scaling functions, 

9iix,Q') = (^)Gi(z^,Q') ^ giix) , 

92{x,Q') = (-^)'g2(i^,Q2) ^ 52(x) , (14) 

which are non-vanishing. 

If the QCD radiative corrections are neglected, gAx, Q^) is related to the 
polarized quark distributions Aq(x) as shown in Eq. (|ll|). In QCD, the distribu- 
tion can be expressed as the Fourier transform of a quark light-cone correlation, 

Aqix,f,') = ge*^-(P5|V^(0)C/(0,An) flj,^(Xn)\PS) (15) 

where n is a light-cone vector (eg. n = (1, 0, 0, —1)) and fi'^ is a renormalization 
scale. U{0, An) — exp(— n ■ A{fin)d^i) is a path-ordered gauge link making 
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the operator gauge invariant. When QCD radiative corrections are taken into 
account, the relation between gi{x,Q^) and Aq(x,fi^) is more comphcated (see 



Sec. 3.5). When is not too large (< 5 GeV^), one must take into account 
the higher-twist contributions to gi{x, Q"^), which appear as l/Q"^ power correc- 
tions. Some initial theoretical estimates of these power corrections have been 



performed Bl 198 



Integrating the polarized quark distributions over x yields the fraction of the 
nucleon spin carried by quarks, 



AE = / dxY] {Aq^{x) + Aq, 
Jo 



(x)) . (16) 



The individual quark contribution Aq is also called the axial charge because 
it is related to the matrix element of the axial current ^'7^75^' in the nucleon 
state. AS is the singlet axial charge. Because of the axial anomaly, it is a 
scale-dependent quantity. 

1.3.2 The Nucleon Spin Sum Rule 

To understand the spin structure of the nucleon in the framework QCD, we can 
write the QCD angular momentum operator in a gauge-invariant form ||186| 



where 



>^QCD — Jq + Jg , (17) 



Jq = d'^X X XTq 



d^x 



Jg = J d'^x X X {E X B) . (18) 

(The angular momentum operator in a gauge-variant form has also motivated 
a lot of theoretical work, but is unattractive both theoretically and experimen- 



tally 1 253].) The quark and gluon components of the angular momentum are 
generated from the quark and gluon momentum densities Tg and Ex B, respec- 
tively. S is the Dirac spin matrix and the corresponding term is the quark spin 
contribution. D — V — igA is the covariant derivative and the associated term 
is the gauge-invariant quark orbital angular momentum contribution. 

Using the above expression, one can easily construct a sum rule for the spin 
of the nucleon. Consider a nucleon moving in the z direction, and polarized 
in the helicity eigenstate A = 1/2. The total helicity can be evaluated as an 
expectation value of in the nucleon state, 

i = iAI](A*) + L,(^) + J,(A*), (19) 
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where the three terms denote the matrix elements of three parts of the angular 



momentum operator in Eq. 18, The physical significance of each term is obvious, 
modulo the momentum transfer scale and scheme dependence (see Sec. 3.5) 



indicated by /i. There have been attempts to remove the scale dependence in 
AS by subtracting a gluon contribution |Q . Unfortunately, such a subtraction 
is by no means unique. Here we adopt the standard definition of AE(/x) as 
the matrix element of the multiplicatively renormalized quark spin operator. 
As has been discussed above, AS(/i) can be measured from polarized deep- 
inelastic scattering and the measurement of the other terms will be discussed 
in later sections. Note that the individual terms in the above equation are 
independent of the nucleon velocity |18S]. In particular, the equation applies 



when the nucleon is traveling with the speed of light (the infinite momentum 
frame) . 

The scale dependence of the quark and gluon contributions can be calculated 
in perturbativ e QCD. B y studying renormalization of the nonlocal operators, 

one can show mi 



^ ( \ _ ^sjy) 1 / -16 3nF A ( Jq{^J) 



As /i — > oo, there exists a fixed-point solution 

ISrif 



(20) 



Jg(oo) 



216 + 3nf ' 



1 16 

J (oo) = . (21) 

Thus as the nucleon is probed at an infinitely small distance scale, approximately 
one-half of the spin is carried by gluons. A similar result has been obtained 
by Gross and Wilczek in 1974 for the quark and gluon contributions to the 



momentum of the nucleon |164|. Strictly speaking, these results reveal little 
about the nonperturbative structure of bound states. However, experimentally 
it is found that about half of the nucl eon momentum is carried by gluons even 



at relatively low energy scales (see eg. |214| ). Thus the gluon degrees of freedom 



not only play a key role in perturbative QCD, but also are a major component of 
nonperturbative states as expected. An interesting question is then, how much 
of the nucleon spin is carried by the gluons at low energy scales? A solid answer 
from the fundamental theory is not yet available. Balitsky and Ji have made an 
estimate using the QCD sum rule approach [pTl: 



I e 



< uaGu >< uu > 



M,-lGeY)^- (22) 

which yields approximately 0.25. Based on this calculation, the spin structure 
of the nucleon would look approximately like 

- = 0.10(from -AY.) + 0.15(from Lg) + 0.25(from Jg) . (23) 
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Lattice [226 and quark model Q calculations of Jq have yielded similar results. 

While AS has a simple parton interpretation, the gauge-invariant orbital 
angular momentum clearly does not. Since we are addressing the structure of 
the nucleon, it is not required that a physical quantity have a simple parton 
interpretation. The nucleon mass, magnetic moment, and charge radius do not 
have simple parton model explanations. The quark orbital angular momentum 
is related to the transverse momentum of the partons in the nucleon. It is well 
known that transverse momentum effects are beyond the naive parton picture. 
As will be discussed later, however, the orbital angular momentum does have a 
more subtle parton interpretation (see Sec. 0). 

In the literature, there are suggestions that r x (— «V) be considered the 
orbital angular momentum [253|. This quantity is clearly not gau ge i nvariant 
and — «V does not correspond to the velocity in classical mechanics [143|. Under 
scale evolution, this operator mixes with an infinite number of other operators 
in light-cone gauge [ 175 1 . More importantly, there is no known way to measure 
such "orbital angular momentum." 



1.3.3 Gluon Helicity Distribution AG'(x, Q^) 

In a longitudinally-polarized nucleon the polarized gluon distribution AG(a;, Q^) 
contributes to spin-dependent scattering processes and hence various experimen- 
tal spin asymmetries. In QCD, using the infrared factorization of hard process, 
IS.G{x,Q'^) with — 1 < a; < 1 can be expressed as 

AG{x, m') = ^ / ^e'^-(F5|^+"(0)L/(0, \n)F+^{\n)\PS) , (24) 

where Fap — {l/2)eaPfiuF^'' . Because of the charge conjugation property of the 
operator, the gluon distribution is symmetric in x: AG'(x, Q^) — AG(— x, Q^). 

The even moments of /^G{x,Q^) are directly related to matrix elements of 
charge conjugation even local operators. Defining 

j dxx"-^AG{x,Q^):^a„{n^) , (n = 3,5---) (25) 

we have 

(PS\Ft'^°'iDt'^ ■■■iD^'"-HF^'l\PS) = 2a„ (^2)5-^1^^2 ... (26) 

(n = 3,5- • •) . 

The renormalization scale dependence is directly connected to renormalization 
of the local operators. Because Eq. ( p^ involves directly the time variable, it is 
difficult to evaluate the distribution on a lattice. However, the matrix elements 
of local operators are routinely calculated in lattice QCD, hence the moments 
of AG{x, Q^) are, in principle, calculable. 

From the above equations, it is clear that the first-moment {n ~ 1) of AG(a;) 
does not correspond to a gauge-invariant local operator. In the axial gauge 
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n ■ A = 0, the first moment of the nonlocal operator can be reduced to a local 
one, E X A, which can be interpreted as the gluon spin density operator. As a 
resuh, the first moment of AG'(a;,/i^) represents the gluon spin contribution to 
the nucleon spin in the axial gauge. In any other gauge, however, it cannot be 
interpreted as such. Thus one can formally write Jg = AG + Lg in the axial 
gauge, where Lg is then the gluon orbital contribution the nucleon spin. There 
is no known way to measure Lg directly from experiment other than defining it 
as the difference between Jg and AG. 

2 Experimental Overview 

A wide variety of experimental approaches have been applied to the investiga- 
tion of the nucleon's spin structure. The experiments complement each other 
in their kinematic coverage and in their sensitivity to possible systematic errors 
associated with the measured quantities. A summary of the spin structure mea- 
surements is shown in Table ^ where the beams, targets, and typical energies 
are listed for each experiment. The kinematic coverage of each experiment is 
indicated in the table by its average four-momentum transfer (Q^) and Bjorken 
X range (for > 1 GeV^). Also given are the average or typical beam and 
target polarizations as quoted by each experimental group in their respective 
publications (or in their proposals for the experiments that are underway). The 
column labeled / lists the dilution factor, which is the fraction of scattered 
events that result from the polarized atoms of interest, and the column labeled 
£ is an estimate of the total nucleon luminosity (# of nucleons/cm^ x of beam 
particles/s) in units of 10^^ nucleons/cm^/s for each experiment. 

In an effort to eliminate possible sources of unknown systematic error in 
the measurements the experiments have been performed with significantly dif- 
ferent experimental techniques. Examples of the large range of experimental 
parameters for the measurements include variations in the beam polarization 
of 40 — 80%, in the target polarization of 30 — 90% and in the correction for 
dilution of the experimental asymmetry due to unpolarized material of 0.1 — 1. 

We now present an overview of the individual experimental techniques with 
an emphasis on the different approaches taken by the various experiments. 
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Table 1: Summary of High Energy Spin Structure Function Measurements. 



Lab 


Exp. 


Year 


Beam 




X 




Pb 


Target 


Pt 


./ 


C X 10"'^^ 










GeV^ 














cm -s 


SLAC 


E80 


75 


10-16 GeV e" 


2 


0.1 - 


0.5 


85% 


H-butanol 


50% 


0.13 


400 




E130 


80 


16-23 GeV e 


5 


0.1 - 


0.6 


81% 


H-biitanol 


58% 


0.15 


400 




E142 


92 


19-26 GeV e~ 


2 


0.03 - 


0.6 


39% 


3 Ho 


35% 


0.35 


2000 




E143 


93 


10-29 GeV e~ 


3 


0.03 - 


0.8 


85% 


NH3 


70% 


0.15 


1000 


















ND.3 


25% 


0.24 


1000 




E154 


95 


48 GeV e~ 


5 


0.01 - 


0.7 


82%) 


■'He 


38% 


0.55 


3000 




E155 


97 


48 GeV e~ 


5 


0.01 - 


0.9 


81% 


NH3 


90% 


0.15 


1000 


















T ;t~\ 
LlD 


000/ 
22% 


IJ.OO 


1 AAA 
iUUU 




Jilioo 


yy 


oU IjreV e 


Q 
O 


n no 
U.Uz - 


u.y 


o6/o 


IN Ms 


7570 


U.lO 


iUUU 


















T iF» 


zz /o 


U.oO 


IUUU 


CERN 


EMC 


85 


100-200 GeV /x+ 


11 


0.01 


0.7 


79% 


NH3 


78% 


0.16 


0.3 




SMC 


92 


100 GeV /i+ 


4.6 


0.006 - 


- 0.6 


82% 


D-butanol 


35% 


0.19 


0.3 






93 


190 GeV n+ 


10 


0.003 - 


-0.7 


80% 


H-butanol 


86% 


0.12 


0.6 






94-95 










81% 


D-butanol 


50% 


0.20 


0.6 






96 










77% 


NHa 


89% 


0.16 


0.6 


DESY 


HERMES 


95 


28 GeV e+ 


2.5 


0.02 - 


0.6 


55% 


^He 


46% 


1.0 


1 






96-97 










55% 


H 


88% 


1.0 


0.1 






98 


28 GeV e- 








55% 


D 


85% 


1.0 


0.2 






99-00 


28 GeV e+ 








55% 


D 


85% 


1.0 


0.2 


CERN 


COMPASS 


01 


190 GeV At+ 


10 


0.005 - 


-0.6 


80% 


NH3 


90% 


0.16 


3 


















LiD 


40% 


0.50 


3 


BNL 


RHIC 


02 


200 GeV p - p 


~ 100 


0.05 - 


0.6 


70% 


Collider 


70% 


1.0 


2 


DESY 


ZEUS/Hl 


?? 


28 X 800 GeV e - p 


22 


0.00006 


- 0.6 


70% 


Collider 


70% 


1.0 


0.2 



2.1 SLAC Experiments 

The SLAC program has focused on high statistics measurements of the inclusive 
asymmetry. The first pioneering experiments on the proton spin structure were 
performed at SLAC in experiments E80 ||3^ and E130 ||6^. These experiments 
are typical of the experimental approach of the SLAC spin program. Polarized 
electrons are injected into the SLAC linac, accelerated to the full beam energy 
and impinge on fixed targets in End Station A. The polarization of the elec- 
trons is measured at low energies at the injector using Mott scattering and at 
high energies in the End Station using MoUer scattering. Target polarization is 
typically measured using NMR techniques. The scattered electrons are detected 
with magnetic spectrometers where electron identification is usually done with 
Cerenkov detectors and Pb-Glass calorimeters. 

For E80 and E130, electrons were produced by photoionization of ^Li pro- 
duced in an atomic beam source. Electron polarization is produced by Stern- 
Gerlach separation in an inhomogeneous magnetic field. Polarized protons were 
produced by dynamic polarization of solid-state butanol doped with a param- 
agnetic substance. Depolarization effects in the target limited the average beam 
currents to ~ 10 nA. In these experiments a considerable amount of unpolarized 
material is present in the target resulting in a dilution of the physics asymmetry. 
For E80 and E130 this dilution reduced the asymmetry by a factor of ~ 0.15. 

Over the last ten years a second generation of high precision measurements 
have been performed at SLAC. Information on the neutron spin structure has 
been obtained using polarized ^He in experiments E142 and E154 Q. 

Here the polarized '^He behaves approximately as a polarized neutron due to the 
almost complete pairing off of the proton spins. The nuclear correction to the 
neutron asymmetry is estimated to be ~ 5 - 10 %. Beam currents were typically 
.5-2 jiK and the polarization was significantly improved for the E154 experi- 
ment using new developments in strained gallium-arsenide photocathodes 
A schematic diagram of the spectrometers used for E142 is shown in Fig. 

Additional data on the neutron and more precise data on the proton has 
come from E143 |, H and E155 ||, || where both and H polarized 
targets using polarized ammonia (NH3 and ND3) and ^LiD were employed. 
The main difference between these two experiments was again an increase in 
beam energy from 26 - 48 GeV and an increase in polarization from 40 % to 80 
%. 

2.2 CERN Experiments 

Following the early measurements at SLAC, the EMC (European Muon Col- 

, |5^ performed the first measurements at x < 0.1. 
Polarized muon beams were produced by pion decay yielding beam intensities 
of lO'^/i/s. The small energy loss rate of the muons allowed the use of very thick 
targets 1 m) of butanol and methane. The spin structure measurements 
by EMC came at the end of a series of measurements of unpolarized nucleon 
and nuclear structure functions, but the impact of the EMC spin measurements 
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Figure 2: Schematic diagram of SLAC E142/E143 spectrometers 



was significant. Their low x measurements, accessible due to the high energy 
of the muons, suggested the breakdown of the naive parton picture that quarks 
provide essentially all of the spin of the nucleon. 

The SMC (Spin Muon Collaboration) experiment |l|, |l|, |l[ [B) began 
as a dedicated follow-on experiment to the EMC spin measurements using an 
upgraded apparatus. An extensive program of measurements with polarized 
and targets was undertaken over a period of ten years. Improvements in 
target and beam performance provided high precision data on inclusive spin- 
dependent structure functions. The large acceptance of the SMC spectrometer 
in the forward direction (see Fig. ||) allowed them to present the first mea- 
surements of spin structure using semi-inclusive hadron production. As with 
EMC, the high energy of the muon beam provided access to the low x regime 
{x < 0.01). 

A new experiment is underway at CERN whose goal is to provide direct infor- 
mation on the gluon polarization. The COMPASS [|lll[ (COmmon Muon Pro- 
ton Apparatus for Structure and Spectroscopy) experiment will use a large ac- 
ceptance spectrometer with full particle identification to generate a high statis- 
tics sample of charmed particles. Using targets similar to those used in SMC 
and an intense muon beam 10^/i/s) improved measurements of other semi- 
inclusive asymmetries will also be possible. 
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Figure 3: Schematic diagram of SMC spectrometer 



2.3 DESY Experiments 

Using very thin gaseous targets of pure atoms (^H, ^H, ■^He) and very high cur- 
rents (~ 40 mA) of stored, circulating positrons or electrons HERMES (HERa 
MEasurement of Spin) has been taking data at DESY since 1995. HERMES 
is a fixed target experiment that uses the stored beam of the HERA col- 
lider .^The polarization of the beam is achieved through the Sokolov-Ternov ef- 



fect 1 266 1, whereby the beam becomes transversely polarized due to a small-spin 
dependence in the synchrotron radiation emission. The transverse polarization 
is rotated to the longitudinal direction by a spin rotator - a sequence of horizon- 
tal and vertical bending magnets that takes advantage of the g — 2 precession 
of the . The beam polarization is measured with Compton polarimeters |]6^ . 

HERMES has focused its efforts on measurements of semi- inclusive asymme- 
tries, where the scattered is detected in coincidence with a forward hadron. 
This was achieved with a large acceptance magnetic spectrometer |ll[ as shown 
in Fig. ^. Initial measurements allowed some limited pion identification with a 
gas threshold Cerenkov detector and a Pb-glass calorimeter. Since 1998 a Ring 
Imaging Cerenkov (RICH) detector has been in operation allowing full hadron 
identification over most of the momentum acceptance of the spectrometer. 

Up to the present, HERMES has taken data only with a longitudinally po- 
larized target. Future runs will focus on high statistics measurements with a 
transversely polarized target to access eg. transversity (see Sect. 6.1) and 92 
(see Sect. 6.2). 
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Figure 4: Schematic diagram of HERMES spectrometer 



Promising future spin physics options also exist at DESY if polarized protons 
can be injected and accelerated in the HERA ring. The HERA-N |205[ | program 
would use the stored 820 GeV proton beam and a fixed target of gaseous polar- 
ized nucleons. This would allow measurements of quark and gluon polarizations 
at ^/s ~ 50 GeV, complimenting the higher energy measurements possible in 
the RHIC spin program. 

A stored polarized proton beam in HERA would also allow e — p collider 
measurements [12C| with the existing HI and ZEUS detectors. Inclusive po- 
larized DIS could be measured to much higher and lower x than existing 
measurements. This would allow improved extraction of the gluon polarization 
via the scaling violations of the spin-dependent cross section. Heavy quark and 
jet production as well as charged-current vs. neutral-current scattering would 
also allow improved measurements of both quark and gluon polarizations. 

2.4 RHIC Spin Program 

The Relativistic Heavy-Ion Collider (RHIC) |254| ] at the Brookhaven National 
Laboratory recently began operations. This collider was designed to produce 
high luminosity collisions of high-energy heavy ions as a means to search for 
a new state of matter known at the quark-gluon plasma. The design of the 
accelerator also allows the acceleration and collision of high energy beams of 
polarized protons and a fraction of accelerator operations will be devoted to 
spin physics with colliding p — p. Beam polarizations of 70% and center-of-mass 
energies of -y/s = 50 — 500 are expected. 



Two large collider detectors, PHEN IX ||235| and STAR ||l67[, along with 



several smaller experiments, BRAHMS [E7§, PHOBOS E73^ and PP2PP, will 
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participate in the RHIC spin program. As an example a schematic diagram of 
the STAR detector is shown in Fig. ||. Longitudinal beam polarization will be 
available for the PHENIX and STAR detectors enabling measurements of quark 
and gluon spin distributions (see Sects. 4.2 and 5.74). 




3 Total Quark Helicity Distribution 

A large body of data has been accumulated over the past ten years on inclu- 
sive polarized lepton scattering from polarized targets. These data allow the 
extraction of the spin structure functions g\'^{x,Q^) and the nearly model- 
independent determination of the total quark contribution to the nucleon spin 
AS — (Aw -I- Au) -I- (Ad -I- Ad) -I- (As -I- As). Inclusive data combined with as- 
sumptions about flavor symmetry, S't/(3) /, and results from beta decay provide 
some model-dependent information on the individual flavor contributions to the 
nucleon spin. Studies of the dependence of gi(x, Q^) allow a first estimate of 
the gluon spin contribution albeit with fairly large uncertainties. These results 
are discussed in the following sections. 
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3.1 Virtual Photon Asymmetries 

Virtual photon asymmetries can be defined in terms of a heHcity decomposition 
of the virtual photon-nucleon scattering cross sections. For a transversely polar- 
ized virtual photon (eg. with helicity ±1) incident on a longitudinally polarized 
nuclcon there are two helicity cross sections cri and era and the longitudinal 
asymmetry is given by 

(71 — (73 

Ai = 21 . (27) 

2 2 

A-z is a virtual photon asymmetry that results from an interference between 
transverse and longitudinal virtual photon-nucleon amplitudes: 

A, = . (28) 



These virtual photon asymmetries, in general a function of x and Q^, are related 
to the nucleon spin structure functions gi{x, Q^) and ,92(2;, Q^) via 

A , n2-. _ 9i{x,Q^)-j'^g2{x,Q'^) 

AM) = ^-^^ , (29) 

where 7 = 2Mx/^/^. 

These virtual photon asymmetries can be related to measured lepton asym- 
metries through polarization and kinematic factors. The experimental longitu- 
dinal and transverse lepton asymmetries are defined as 



All 



A^ = 11^^, (30) 

where a'^'^ (crTi) jg ^jjg cross section for the lepton and nucleon spins aligned 
(anti-aligned) longitudinally, while cr^^ (^T^) the cross section for longitudi- 
nally polarized lepton and transversely polarized nucleon. The lepton asymme- 
tries are then given in terms of the virtual photon asymmetries through 

A|| = D{Ai+r]A2) 

Ai. = rf(A2-CAi) . (31) 

The virtual photon (de)polarization factor D is approximately equal to y = u/E 
(where u is the energy of the virtual photon and E is the lepton energy), but is 
given explicitly as 

D=[l-{l-y)e]/{l + eR) , (32) 
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where e is the magnitude of the virtual photon's transverse polarization 

e = [4(1 -y)- i^i']/[2y^ + 4(1 - y) + ^^y"] . (33) 

and 

R = (tl/^t (34) 

is the ratio of longitudinal to transverse virtual photon cross sections. 
The other factors are given by 

77 = e72//[l - e(l - v)] (35) 



^-^Vlf^ (36) 
3.2 Extraction of Q^) 

The nucleon structure function is extracted from measurements of the lepton- 
nucleon longitudinal asymmetry (with longitudinally polarized beam and target) 

^11 = ^Ti+aTT ' (38) 

where a^^ {a^'^) represents the cross section when the electron and nucleon spins 
are aligned (anti-aligned). These cross sections can also be expressed in terms 
of spin- independent ajj and spin-dependent up cross sections 



(T 



TT 



(Ju + <^p (39) 
cr^i ^ au -(TP . (40) 

In the limit of stable beam currents, target densities and polarizations, the 
experimentally measured asymmetry Aoxp is usually expressed in terms of the 
measured count rates N and the number of incident electrons Nb 



^11 is then determined via 



where Pp and Px are the beam and target polarizations respectively, / is a 
dilution factor due to scattering from unpolarized material and Apc accounts 
for QED radiative effects |1. 
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If however there is a time variation of the beam or target polarization or 
luminosity, the asymmetry should be determined using 

NUlII _ AtTTt^TX 
^11 = fT fi- + ^RC , (43) 

since in this case the measured count rates can be written in terms of au and 
crp 

TVTT = auf n^J{t)dt + apf n^^ {t)PB{t)PT{t)dt = auL^^ + <JpL^^ 

ATTi ^ auJnl^{t)dt-apJn\^{t)PB{t)PT{t)dt = auL^^+apL^p^, ^ ' 

where now ub represents the product of beam current and target areal density - 
the luminosity. In Eq. ^ we have ignored a factor accounting for the acceptance 
and solid angle of the apparatus which is assumed to be independent of time. 

The spin structure function gi{x, Q^) can then be determined from the lon- 
gitudinal asymmetry A||(a;, Q^), 

ffi= + (7-^7)^2] , (45) 

where Fi = Fi(x,Q^) is the unpolarized structure function. The unpolarized 
structure function Fi is usually determined from measurements of the unpolar- 
ized structure function F2 and R using 

-Fi -i^2(l+7')/(2a:(l + i?)) . (46) 
To use the above equation we need an estimate for 1 is constrained 



to be less th an V-R |124], but A2 can also be determined from measurements 
(see Sec. 6.1) with a longitudinally polarized lepton beam and a transversely 
polarized nucleon target (when combined with the longitudinal asymmetry). 

As a guide to the relative importance of various kinematic terms in the above 
equations we present examples of the magnitude of these terms in Table ||, 
typical for the SMC and HERMES experiments. 

For extraction of the neutron structure function gi from nuclear targets, 
eg. and '^He, additional corrections must be applied. For the deuteron, the 
largest contribution is due to the polarized proton in the polarized deuteron 
which must be subtracted. In addition a D-state admixture into the p — n wave 
function will reduce the deuteron spin structure function due to the opposite 
alignment of the p — n spin system in this orbital state; thus 



^2, 2.gf(x,Q2) 



9ii^.Q') = (1 i.sl,^) -gi(^.Q^) (47) 
where ujo is the D-state probability of the deuteron. Typically a value of loo = 



0.05 ± 0.01 1 211 is used for this correction. 

For polarized "^He, a wavefunction correction for the neutron and proton 
polarizations is applied using 

gnx,Q') = -{gl"'-2p,gP,) , (48) 
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Tabic 2: Typical kinematic factors entering into the extraction of gi. Examples 
from the SMC and HERMES experiments are given. 



SMC 




< a; > 


< 0" > 


y 


7 


e 


7 


- V 


D 


0.005 


1.30 





729 





008 





505 





005 





721 


0.008 


2.10 





736 





010 





493 





006 





745 


0.014 


3.60 





721 





014 





517 





008 





748 


0.025 


5.70 





639 





020 





638 





009 





671 


0.035 


7.80 





625 





024 





657 





Oil 





666 


0.049 


10.40 





595 





029 





695 





012 





643 


0.077 


14.90 





543 





037 





756 





014 





592 


0.122 


21.30 





490 





050 





809 





016 





545 


0.173 


27.80 





451 





062 





843 





018 





508 


0.242 


35.60 





413 





076 





873 





020 





468 


0.342 


45.90 





376 





095 





897 





022 





428 


0.480 


58.00 





339 





118 





919 





024 





384 


HERMES 




<x> 


<Q'> 


y 


7 


e 


7 


- V 


D 


0.023 


0.92 





775 





045 





427 





029 





778 


0.033 


1.11 





652 





059 





620 





028 





635 


0.047 


1.39 





573 





075 





721 





030 





547 


0.067 


1.73 





500 





096 





798 





032 





476 


0.095 


2.09 





426 





123 





861 





034 





405 


0.136 


2.44 





348 





163 





913 





035 





329 


0.193 


2.81 





282 





216 





945 





037 





268 


0.274 


3.35 





237 





281 





962 





040 





227 


0.389 


4.25 





212 





354 





969 





047 





208 


0.464 


4.80 





200 





397 





972 





050 





198 


0.550 


5.51 





194 





440 





973 





055 





195 


0.660 


7.36 





216 





457 





965 





065 





224 



20 



where p„ = (86 ± 2)% and pp = (—2.8 ± 0.4)% as taken from a number of 
calculations [ 146| , 106{ . Additional corrections due to the neutron binding energy 
and Fermi motion have also been investigated 
relatively small. 



a, 106: 26C| and shown to be 



3.3 Recent Results for gi{x,Q'^) 

Most of the experiments listed in Table |^ have contributed high precision data 
on the spin structure function gi^XjQ"^). Where there is overlap (in x and Q^), 
the agreement between the experiments is extremely good. This can be seen in 
Fig. pi where the ratio of the polarized to unpolarized proton structure function 
g\/F[ is shown. Analysis of the dependence of this ratio [Q has shown that 
it is consistent experimentally with being independent of within the range 
of existing experiments, although this behavior is not expected to persist for all 
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Figure 6: Ratio of polarized to unpolarized proton structure function from the 
SMC, E143 and HERMES experiments. 

A comparison of the spin structure functions gP'"^'" are shown in Fig. |^. Some 
residual dependence is visible in the comparison of the SMC data with the 
other experiments. The general dependence of gi will be discussed in Sect. 
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3.5. 
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Figure 7: Compilation of the world's data on g^^gf^Qi 

While the results shown in Figs. correspond to > 1 GeV^, data also 
exists at lower Q^, because of the large kinematic acceptance in many of the 



experiments. Much of this data 0, Ea, 238 1, when expressed as gi/Fi appears 
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to be largely independent of Q^. 



3.4 First Moments of gi{x,Q'^) 

The initial interest in measurements of gi{x, Q^) was in comparing the measure- 
ments to several predicted sum rules, specifically the Ellis- Jaffe and Bjorken sum 
rules. These sum rules relate integrals over the measured structure functions to 
measurements of neutron and hyperon beta-decay. 



The Ellis- Jaffe sum rule [137| starts with the leading-order QPM result for 
the integral of 51(2;): 



/ gi{x)dx = l [ y2e^Aq,{x)dx , (49) 
Jo ^ Jo ^ 

where the sum is over u,d,s,u,d,s for three active quark flavors and the 
dependence has been suppressed as it is absent in the simple QPM. Introducing 
the SU{3)f nucleon axial charges: 

ao = {Au + As) + (Ad -I- Ad) + (As -I- As) (50) 
03 = (Au + Au) - (Ad + Ad) (51) 

as = (Am + Au) + (Ad + Ad) - 2(As + As) , (52) 

where Aq^ = Aqi{x)dx, the Ellis- Jaffe sum rule then assumes that the strange 
quark and sea polarizations are zero (As = A^j = 0). Then for the proton and 
neutron integrals the Ellis- Jaffe sum rule gives: 

ri= f g'l{x)= ^03 + 4^8 + 4ao = 0.186 ±0.004 









^"3 - 


f ^ag - 


f ^flo 












-^03 4 







= [ gP{x) = -^03 + ^ag -I- ^ao = -0.025 ± 0.004 . (53) 
"'0 

To evaluate the integrals it is assumed that oq = ag which is true if As = 0. 
Then 03 is determined from the ratio of axial- vector to vector coupling constants 
in neutron decay = —gA/gv = 1.2670 ± 0.0035 |24l| . A value for ag can 
be estimated with the additional assumption of SU(3) flavor symmetry which 
allows one to express gA/gv for hyperon beta decays in terms of 03 and ag (see 
Table ||), giving ag — 0.58 ±0.03. Nucleon and hyperon beta decay is sometimes 
parameterized in terms of the F and D coefficients. These coefflcients are related 
to the axial charges 03 and ag with 

as = F + D 

ag 3F-D . (54) 

The assumptions implicit in the Ellis- Jaffe sum rule, eg. As = Aqi = and 
SU (3) / symmetry, may be signiflcantly violated. On the contrary, the Bjorken 
sum rule 

1 

/ [gf{x)-g'^{x)]dx^ -as = 0.211 ±0.001 (55) 
Jo o 
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Table 3: Relation of the neutron and hyperon beta decays to the nucleon's axial 
charges (as defined in the text) assuming SU{3)f symmetry. 



Decay 



qaI gv in terms of Axial charges 



Experimental |241 



n - 


-> pe_ 


Ve 


A- 




Ve 




— > ne 











-as 



-1.2670 ±0.0035 
-0.718 ±0.015 
0.340 ±0.017 
-0.25 ±0.05 



requires only current algebra and isospin symmetry (eg. Am^ = Ad") in its 
derivation. Note that both the Ellis- Jaffe and Bjorken sum rules must be cor- 
rected for QCD radiative corrections. For example, these corrections have been 
evaluated up to order [217| and amount to ~ 10% correction for the EUis- 
Jaffe sum rule and ~ 15% correction for the Bjorken sum rule at = 5 GeV^. 

Comparison of these predictions with experiment requires forming the inte- 
grals of the measurements of 171(3;, Q^) over the full x range from ^ 1 at a fixed 

. Thus extrapolations are necessary in order to include regions of unmeasured 
X, both at high and low x. For the large x region this is straightforward: since 
gi{x) is proportional to a difference of quark distributions it must approach 
zero as X —> 1 as this is the observed behavior of the unpolarized distribu- 
tions. However the low x region is problematic, as there is no clear dependence 
expected. In the first analyses simple extrapolations based on Regge parame- 
terizations |17C, 135| were used. Thus gi{x) was assumed to be nearly constant 
for X — > 0. Later, Next-to-Leading-Order (NLO) QCD calculations |136] (see 
Sec. 3.5) suggested that these parameterizations likely underestimated the low 
X contributions. The NLO calculations cannot predict the actual x dependence 
of the structure function, but can only take a given x dependence and predict its 
dependence on . Thus by using the Regge parameterizations for low ^ 1, 
they can give the low x behavior at the of the experiments, eg. 3 — 10 

GeV2. 

Evaluating the experimental integrals at a fixed requires an extrapola- 
tion of the measured structure function. In general, for each experiment, the 
experimental acceptance imposes a correlation between x and preventing a 
single experiment from measuring the full range in a; at a constant value of . 
Thus the data must be QCD-evolved to a fixed value of . This has often 
been done by exploiting the observed independence of gi{x, Q'^)/Fi{x, Q"^) 
(see Fig. ||). In this case most of the dependence of gi{x,Q^) results from 
the dependence of the unpolarized structure function Fi{x, Q^) which is well 
measured in other experiments. Alternatively, NLO QCD fits (as described in 
the next section) can be used to evolve the data sets to a common Q^. 

The El 55 collaboration has recently reported |^ a global analysis of spin 
structure function integrals. They have evolved the world data set on (/f (x, Q^) 
and gi{x,Q^) to = 5 GeV^ and have extrapolated to low and to high x 
using a NLO fit to the data. Their results are compared in Table ^ with the 



24 



predictions for the Ellis- Jaffe and Bjorken sum rules (Eqs. ^3|p5| ) including QCD 
radiative corrections for — h GeV^ up to order using the calculations of 
Ref. 1 217 and world-average for [241|. 



Table 4: Comparison of Sum Rule predictions including corrections up to order 



with a global analysis of the experiments. 



Sum Rule 


Calculation 


Experiment |49|] 


EJ Sum rf (g^ = 5 GeV^) 
EJ Sum r'/(Q2 = 5 GeV^) 
Bj rP-r?(Q2 = 5 GeV2) 


0.163 ±0.004 
-0.019 ±0.004 
0.181 ±0.005 


0.118 ±0.004 ±0.007 
-0.058 ± 0.005 ±0.008 
0.176 ±0.003 ±0.007 



As seen in Table | the Bjorken sum rule is well verified. In fact some anal- 



yses ||136[ have assumed the validity of the Bjorken sum rule and used the 
dependence of T\ — F" to extract a useful value for a^. In contrast there is a 
strong violation of the Ellis- Jaffe sum rules. Many early analyses of these re- 
sults interpreted the violation in terms of a non-zero value for As (in which case 
Oq ^ ag), using only the leading order QPM. However modern analyses have 
demonstrated that a full NLO analysis is necessary in order to interpret the re- 
sults. This analysis will be described in the next section. Here, for completeness, 
we give the leading order QPM result. 

Within the leading order QPM, Au + Au, Ac?+ Ad and As + As can be de- 
termined by using Eqs. 53 with the experimental values from Table ^ Dropping 
the assumption of aq — ag, but retaining the SU{S)f assumption to determine 
ag, one finds: 

Au + Au = 0.78±0.03, Ad+Ad = -0.48±0.03, As + As = -0.14±0.03 (56) 

after applying the relevant QCD radiative corrections to the terms in Eq. ^ 
(corresponding to a factor of 0.859 multiplying the triplet and octet ag charges 
and a factor of 0.878 multiplying the singlet ao charge for = 5 GeV^). This 
then gives a very small value for the total quark contribution to the nucleon's 
spin, AE — 0.16 ± 0.08. Note that the quoted uncertainties reflect only the 
uncertainty in the measured value of Fi and not possible systematic effects due 
to the assumption of SU (3) / symmetry and NLO effects. Studies of the effect 
of SU{3)f symmetry violations have been estimated ^ to have little effect on 
the uncertainty in Au and Ad, but can increase the uncertainty on As by a 
factor of two to three. NLO effects are the subject of the next section. 



3.5 Next-to-Leading Order Evolution of gi{x,Q^) 

As discussed above the spin structure functions possess a significant depen- 
dence due to QCD radiative effects. It is important to understand these effects 
for a number of reasons, including comparison of different experiments, forming 
structure function integrals, parameterizing the data and obtaining sensitivity 
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to the gluon spin distribution. As the experiments are taken at different ac- 
celerator facilities with differing beam energies the data span a range of . 
In addition, because of the extensive data set that has been accumulated and 
the recently computed higher-order QCD corrections, it is possible to produce 
parameterizations of the data based on Next-to-Leading-Order (NLO) QCD fits 
to the data. This provides important input to future experiments utilizing po- 
larized beams (eg. the RHIC spin program). These fits have also yielded some 
initial information on the gluon spin distribution, because of the radiative effects 
that couple the quark and gluon spin distributions at NLO. 

At NLO the QPM expression for the spin structure function becomes 

gx{x,Q^) = \^elCq{x,as)® t\q^x,Q^) + ^Cg{x,as)® t^G{x,Q^) , (57) 

where for three active quark flavors {Nf — 3) the sum is again over quarks 
and antiquarks: u,d, s,u,d,s. Cq{x,as) and Cg(x, a^) are Wilson coefficients 
and correspond to the polarized photon-quark and photon-gluon hard scattering 
cross section respectively. The convolution (g) is defined as 

C{x, as) ® q{x, Q2) = f ^c{-,as)q{x, Q^) . (58) 

Jx y V 

The explicit dependence of the nucleon spin structure function on the gluon spin 
distribution is apparent in Eq. At Leading Order (LO) C° — (5(1 — x) and 
Cg = and the usual dependence (Eq. ^9|) of the spin structure function on the 
quark spin distributions emerges. At NLO however, the factorization between 
the quark spin distributions and coefhcient functions shown in Eq. |5^ cannot be 
defined unambiguously. This is known as factorization scheme dependence and 
results from an ambiguity in how the perturbative physics is divided between 
the definition of the quark/gluon spin distributions and the coefficient functions. 
There arc also ambiguities associated with the definition of the 75 matrix in n 



dimensions | 272 | and in how to include the axial anomaly. This has lead to a 
variety of factorization schemes that deal with these ambiguities by different 
means. 

We can classify the factorization schemes in terms of their treatment of 
the higher order terms in the expansion of the coefficient functions. The 
dependence of this expansion can be written as: 

Q(x,a,) = C°(x)-f^^C«-t-... . (59) 



In the so-called Modified-Minimal-Subtraction (MS) scheme |231, 277 1 the first 
moment of the NLO correction to Cg vanishes (i.e. jj^ c\^\x)dx = 0), such that 
AG does not contribute to the first moment of gi . In the Adler-Bardeen ||6^, p8| 
scheme (AB) the treatment of the axial anomaly causes the first moment of Cg 
to be non-zero, leading to a dependence of J gi{x)dx on J AG{x)dx. This then 



26 



leads to a difference in the singlet quark distribution AE in the two schemes: 



^'^ J X y 

AGix,Q^)AB = AG(:E,g2) (60) 



A third scheme, s omet imes called the JET scheme 101 , 218 or chirally in- 



variant (CI) scheme 104 , is also used. This scheme attempts to include all 
perturbative anomaly effects into Cg. Of course any physical observables (eg. 
gi{x,Q'^)) are independen t of t he c hoice of scheme. There are also straightfor- 
ward transformations ^37[ ^19| that relate the schemes and their results to 
one another. 

Once a choice of scheme is made the dependence of gi can be calculated 



using the Dokahitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [163| equations. 
These equations characterize the evolution of the spin distributions in terms of 
Q^-dependent splitting functions Pij{x, as): 



AqNsix,Q') = ^^P:^(^^Aq_ 



d f AJ: \ a,(Q2) 



NS 



dlnQ2 V AG y 27r 









\ Pgq 


^99 J 





(61) 



where the non-singlet quark distributions Aqjv5'(a;, Q^) for three quark flavors 
are defined with 

AqNs{x, Q^) = (Au -I- Au) - ^(Ad + Ad) - i(As + As) . (62) 
The splitting functions Pij can be expanded in a form similar to that for the 



coefficient functions Gi(a;, a^) in Eq. 59 and have been recently evaluated [231 



277| in NLO. 

The remaining ingredients in providing a fit to the data are the choice of 
starting momentum scale Qq and the form of the parton distributions at this 
Qq. The momentum scale is usually chosen to be < 1 GeV^ so that the quark 
spin distributions are dominated by the valence quarks and the gluon spin dis- 
tribution is likely to be small. Also, as discussed above, at lower momentum 
transfer some models for the x dependence of the distributions (eg. Regge-type 
models for the low x region) are more reliable. The form of the polarized parton 
distributions at the starting momentum scale are parameterized by a variety of 
X dependences with various powers. This parameterization is the source of some 
of the largest uncertainties as the x dependence at low values of a; < 0.003 is 
largely unconstrained by the measurements. As an example, Ref. ||3^ assumes 
for one of its fits that the polarized parton distributions can be parameterized 

by 

Aq,{x,Ql)^ A,x''' {l-xf^ {l + -f,x'^) . (63) 

With such a large number of parameters it is usually required to place additional 
constraints on some of the parameters. Often SU (3) / symmetry is used to 
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constrain the parameters, or the positivity of the distributions ( | Ag^ (a;) | < qi{x)) 
is enforced (note that this positivity is strictly vaUd only when all orders are 
included; see Ref. |3^). Thus in other fits, the polarized distributions are taken 
to be proportional to the unpolarized distributions as in eg. Ref. p9[: 



(64) 



A larg e number of NLO f it s ha ve r ecent ly been published |15£, 151 

IS 



16g, gl9|, g20|, g2l|, ^ |49|, |llj. These fits include a wide variety 



of assumptions for the forms of the polarized parton distributions, differences 
in factorization scheme and what data sets they include in the fit (only the 
most recent fits [161| include all the published inclusive data). Some fits [115| 
have even performed a NLO analysis including information from semi-inclusive 
scattering (see Sec. 4.1). A comparison of the results from some of these recent 
fits is shown in Table ^. 
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Table 5: Results from NLO fits to data for first moments of quark and gluon distributions. Missing data refers to data sets 
that are not included in the fits, PPDF refers to the assumptions for Polarized Parton Distribution Functions and Q^^ refers 
to the evolved where the first moments are evaluated. 



Reference 


Scheme 




Missing Data 


PPDF 




AS 


AG 






GeV2 






GeV2 






ABFR98 [38] 


AB 


1 


HERMES (p) 




1 


0.41 ±0.03 


0.95 ±0.18 


(Fit-A) 






E155(pd) 
















Semi-inc 










LSS99 [221] 


JET 


1 


Semi-inc 


Aqi cx qi 


1 


0.39 ±0.04 


0.57±0.14 




AB 


1 








0.41 ±0.04 


0.58 ±0.04 




MS 


1 








0.28 ±0.04 


0.07 ±0.10 


GOTOOO [161] 


MS 


1 


Semi-inc 


Aqi cx qi 


1 


0.050 


0.53 


(NLO-1) 




1 






5 


0.054 


0.86 






1 






10 


0.055 


1.0 


FSOO [115] 


MS 


0.5 




Aqi cx qi 


10 


0.050 


0.53 


(ii) 

















n □ □ □ 



Note that in the JET and AB schemes AS includes a contribution from AG. 
Thus the overriding result of these fits is that the quark spin distribution AS 
is constrained between 0.05 — 0.30 but that the gluon distribution and its first 
moment are largely unconstrained. The extracted value for AG{Q^ = 5GeV^) 
is typically positive but the corresponding uncertainty is often 50 — 100% of 
the value. Note that the uncertainties listed in Table || are dependent on the 
assumptions used in the fits. 

Estimates of the contribution from higher twist effects 198 1 cor- 
rections) suggest that the effects are relatively small at the present experimental 
Q^. This is further supported by the generally good fits that the NLO QCD 
calculations can achieve without including possible higher-twist effects. 

Lattice QCD calculations of the first mom e nts and secon d moments of the 
polarized spin distributions are underway |15C, 125, 158, 165|. Agreement with 
NLO fits to the data is reasonable for the quark contribution, although the 
Lattice calculations are not yet able to calculate the gluon contribution. 



4 Individual Quark Helicity Distributions 

As shown in the last section, the inclusive lepton asymmetries generally provide 
spin structure information only for the sum over quark flavors. Access to the in- 
dividual flavor contributions to the nucleon spin requires assumptions including 
SU (3) / symmetry in the weak decay of the octet baryons (nucleons and strange 
hyperons) . 

Potentially more direct information on the individual contributions of u, d, 
and s quarks as well as the separate contributions of valence and sea quarks is 
possible via semi-inclusive scattering. Here one or more hadrons in coincidence 
with the scattered lepton are detected. The charge of the hadron and its valence 
quark composition provide sensitivity to the flavor of the struck quark within 
the Quark-Parton Model (QPM). 

Semi-inclusive asymmetries also allow access to the third leading-order quark 
distribution 6q called transversity. Because of the chiral odd structure of this 
distribution function it is not measurable in inclusive DIS. Transversity will be 



discussed in Sec. 3.2. Additionally, semi-inclusive asymmetries can provide a 
degree of selectivity for different reaction mechanisms that are sensitive to the 
gluon polarization. The sensitivity of semi-inclusive asymmetries to the gluon 
polarization will be discussed in Sect. ^. The flavor decomposition of the nucleon 
spin using semi-inclusive scattering will be discussed in the next two sections. 

4.1 Semi-Inclusive Polarized Lepton Scattering 

Within the QPM, the cross section for leptoproduction of a hadron (semi- 
inclusive scattering) can be expressed as 



-r- = crDis(x,Q ) 
az 



j:^elqdx,Q^)Dliz,Q^) 



(65) 
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where ctdis is the inclusive DIS cross section, the fragmentation function, 
D^{z, Q^), is the probability that the hadron h originated from the struck quark 
of flavor i, z = Eh/i' is the hadron momentum fraction and the sums are over 
quark and antiquark flavors u,d,s,u,d,s. To maximize the sensitivity to the 
struck "current" quark, kinematic cuts are imposed on the data in order to 
suppress effects from target fragmentation. These cuts typically correspond to 
>9-l0 GeV2 and z > 0.2. 

In general the fragmentation functions D'^{z, Q^) depend on both the quark 
flavor and the hadron type. In particular for a given hadron ^ D^. This 
effect can be understood in terms of the QPM: if the struck quark is a valence 
quark for a particular hadron, it is more likely to fragment into that hadron (eg. 
D'!^^ > ). A flavor sensitivity is therefore obtained as is a sensitivity to the 
antiquarks (eg. ZJj^ > D^^). 

Eq. ^ displays a factorization of the cross section into separate z and x 
dependent terms. This is an assumption of the QPM and must be experimen- 
tally tested. Measurements of unpolarized hadron leptoproduction have 
shown good agreement with the factorization hypothesis. Data from e"*" — e~ — > 



hadrons can also be used to extract fragmentation functions 209 . Both the 
and z dependence of the fragmentat ion functions have been parameterized 



within string models of fragmentation |263| that are in reasonable agreement 
with the measurements. Recently the dependence of the fragmentation func- 
tions have been calculated to NLO |73| . 

Assuming factorization of the cross section as given in Eq. 
the asymmetry for leptoproduction of a hadron as 

e^^Aq,{x,Q^)J^D'^{z,Q^)dz 
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we can write 



J:,eU^{^,Q^)J^D'^iz,Q^)dz 



(66) 



Due to parity conservation the fragmentation functions contain no spin de- 
pendence as long as the final-state polarization of the hadron is not measured 
(spin dependent fr agm entation can be accessed through the self-analyzing de- 
cay of A - see Sec. |8.2| ). By making measurements with H,D and '^He targets 
for different final-state hadrons and assuming isospin symmetry of the quark 
distributions and fragmentation functions a system of linear equations can be 
constructed: 

\ 



AP 

7T 

An 
An 



Ad 

As 
Au 



(67) 



V : / V : / 

and solved for the Aqi. In these equations, the unpolarized quark distributions 
are taken from a variety of parameterizations (eg. Ref. |157| , ^14[ ) and the 
fragmentation functions are taken from measurements |m, |209|] or parameteri- 



zations 1 263 



EMC, SMC and HERMES have made measurements of semi-inclusive asym- 
metries. A comparison of the measurements from SMC and HERMES is shown 
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Figure 8: Virtual photon asymmetries for semi-inclusive lepton scattering. In- 
clusive asymmetries are shown in the leftmost panels for comparison. 



in Fig. |. As the HERMES data are taken at < 2.5 GeV^ and the 

SMC data at < >— 10 GeV^, these data suggest that the semi-inclusive 
asymmetries are also approximately independent of . 

It is important to note, especially for the lower data of HERMES, that 
Eq. ^ must be modified if parameterizations of the unpolarized quark distri- 
butions are used. In some parameterizations it is assumed that the unpolarized 
structure functions are related by the Callen-Gross approximation Fi — F2/2X 
rather than by the complete expression Fi = -^2(1 + 7^)/ (2a;(l Thus some 

experimental groups will present Eq. ^ with an extra factor of (H- -I- 7^) 

included. 

Up to now results have only been reported for positively and negatively 
charged hadrons (summing over tt, p and K) because of the lack of sufficient 
particle identification in the experiments. This reduces the sensitivity to some 
quark flavors (eg. strangeness) and requires additional assumptions about the 
flavor dependence of the sea quark and anti-quark distributions. Two assump- 
tions have been used to extract information on the flavor and sea dependence 
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of the quark polarizations, namely 



Ausix) _ Ads{x) _ As{x) _ Au{x) _ Ad{x) _ As{x) 
Us{x) ds{x) s{x) u{x) d{x) s{x) 

or 

Aus{x) = Ads{x) = As{x) = Au{x) = Ad{x) = As{x) . (69) 

Here Aus and Ads represent the u and d sea quark spin distributions. A 
comparison of the extracted valence and sea quark distributions from HER- 
MES and SMC is shown in Fig. ^. The valence distributions are defined using 
Aqy — Aq — Aq. Typical systematic errors are also shown in Fig. |^ and include 
the difference due to the two assumptions for the sea distributions given by 
Eqs. |6^ The solid hues are positivity limits corresponding to Aq{x) = q{x). 
The da shed lines are parameterizations from Gehrmann and Stirling (Gluon 



A-LO) ll51 



Values for the integrals over the spin distributions from SMC and HERMES 
are compared in Table |6[ The dominant sensitivity to Au within the quark sea 
is due to the factor of two larger charge compared to d and s. 



Table 6: Comparison of the first moment of separated quark spin distributions 
as determined from semi-inclusive DIS lepton scattering. Both statistical and 
systematic uncertainties are given. 



/g Aqi{x)dx 


SMC results 
Q2 = 10 GeV^ 


HERMES results 
Q2 = 2.5 GeV^ 


Auv 
Ady 
Au 


0.77 ±0.10 ±0.08 
-0.52 ±0.14 ±0.09 
0.01 ±0.04 ±0.03 


0.57 ±0.05 ±0.08 
-0.22 ±0.11 ±0.13 
-0.01 ±0.02 ±0.03 



While the experimental results presented in Table ^ have been extra cted 
through a Leading-Order QCD analysis, NLO analyses are possible [114| and 
several such analyses have recently been published |ll5| . 

Future measurements from HERMES and COMPASS will include full par- 
ticle identification providing greater sensitivity to the flavor separation of the 
quark spin distributions. In particular, due to the presence of strange quarks in 
the K valence quark distribution, K identification is expected to give significant 
sensitity to As{x). 



4.2 High Energy p — p Collisions 

The production of weak bosons in high energy p — p collisions at RHIC 
provides unique sensitivity to the quark and antiquark spin distributions. The 
maximal parity violation in the interaction and the dependence of the produc- 
tion on the weak charge of the quarks can be used in principle to select specific 
flavor and charge for the quarks. Thus the single spin longitudinal asymmetry 
for production {pp W^X) can be written |^ 
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• HERMES 

o SMC 




Figure 9: Quark flavor spin structure from hadron leptoproduction from the 
HERMES |l| and SMC ||^, || experiments. Both data sets have been evolved 
to < >— 2.5 GeV^. Typical systematic errors from Ref. are shown 
in the shaded band. The solid lines are the positivity limits corresponding to 
Aq[x) = q{x). The dashed lines are the parameterization from Gehrmann and 



Stirling (Gluon A-LO) [151| 
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^ Au{xi)d{x2) - Ad{xi)u{x2) 
u{xi)d{x2) + d{xi)u{x2) 

where xi and X2 refer to the x value of the quark and antiquark participating in 
the interaction (see for example Fig. [l0| ). Making the replacement u <-> d gives 
the asymmetry for W~ production. In the experiments the are detected 
through their decay to a charged lepton {fj,^ in PHENIX and in STAR) 
and the xi,X2 values are determined from the angles and energies of those 
detected leptons. Thus for production with xi ^ X2 the valence quarks 
are selected for xi and Al{W'^) ^ Au{xi)/u{xi), while for xi ^ X2 valence 
quarks are selected for X2 and Al{W^) ~ Ad{xi) /d{xi). Detection of then 
gives Au/u and Ad/d. An example of the expected sensitivity of the PHENIX 
experiment after about four years of data taking is shown in Fig. 



5 Gluon Helicity Distribution 

As remarked in the Introduction, the gluon contribution to the spin of the nu- 
cleon can be separated into spin and orbital parts. As with its unpolarized 
counterpart, the polarized gluon distribution is difficult to access experimen- 
tally. There exists no theoretically clean and, at the same time, experimentally 
straightforward hard scattering process to directly measure the distribution. In 
the last decade, many interesting ideas have been proposed and some have led 
to useful initial results from the present generation of experiments; others will 
be tested soon at various facilities around the world. 

In the following subsections, we discuss a few representative hard-scattering 
processes in which the gluon spin distribution can be measured. 




Figure 10: W boson production in pp collisions. 
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Figure 11: Predicted sensitivity of the PHENIX detector at RHIC for 
ment of the quark flavor spin contributions. 
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Figure 12: Typical gluon helicity distributions [151| obtained from fits to avail- 
able polarized DIS data. 



5.1 AG{x) from QCD Scale Evolution 

As discussed in Section 3.5, the polarized gluon distribution enters in the factor- 
ization formula for spin-dependent inclusive deep-inelastic scattering. Since the 
51 (cc, Q^) structure function involves both the singlet quark and gluon distribu- 
tions as shown in Eq. only the dependence of the data can be exploited 
to separate them. The dependence results from two different sources: the 
running coupling as{Q'^) in the coefficient functions and the scale evolution of 
the parton distributions. As the gluon contribution has its own characteristic 
behavior, it can be isolated in principle from data taken over a wide range 
of Q2. 

Because the currently available experimental data have rather limited 
coverage, there presently is a large uncertainty in extracting the polarized gluon 
distribution. As described in Sec. 3.5, a number of NLO fits to the world 
data have been performed to extract the polarized parton densities. While 
the results for the polarized quark densities are relatively stable, the extracted 
polarized gluon distribution depends strongly on the assumptions made about 
the x-dependence of the initial parameterization. Different fits produce results 
at a fixed x differing by an order of magnitude and even the sign is not well 
constrained. 

Several sets of polarized gluon distributions have been used widely in the 
literature for the p urpo se of estimating outcomes for future experiments. An 
example from Ref. [ 151 1 of the range of possible distributions is shown in Fig. 
Of course the actual gluon distribution could be very different from any of these. 



5.2 AG{x) from Di-jet Production in e — p Scattering 

In lepton-nucleon deep-inelastic scattering, the virtual photon can produce two 
jets with large transverse momenta from the nucleon target. To leading-order 
in as, the underlying hard scattering subprocesses are Photon-Gluon Fusion 
(PGF) and QCD Compton Scattering (QCDC) as shown in Fig. |l|. If the initial 
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b) 



Figure 13: Leading-order Feynman diagrams for di-jet production in DIS: (a) 
Photon-Gluon Fusion, (b) Photon-Quark Compton scattering. 



photon has momentum q and the parton from the nucleon (with momentum P) 
has momentum xP, the invariant mass of the di-jet \s s = {q + xP)^, the x at 
which the parton densities are probed is 

= -B (^) . (71) 

where is the Bjorken x variable. Therefore the di-jet invariant mass fixes 
the parton momentum fraction. Depending on the relative sizes of s and Q^, 
Xp can be an order of magnitude larger than xb- 

If the contribution from the quark initiated subprocess is small or the quark 
distribution is known, the two-jet production is a useful process to measure the 
gluon distribution. The di-jet invariant mass provides direct control over the 
fraction of the nucleon momentum carried by the gluon {xq — xp). Indeed, 
di-jet data from HERA have been used by the HI and ZEUS collaborations to 
extract the unpolarized gluon distribution ||2^, |255| . With a polarized beam and 
target, the process is ideal for probing the polarized gluon distribution. 

The un pola rized di-jet cross section for photon-nucleon collisions can be 



written as [ |118 



adiHet (<Z, xP) = aPGF , + = AG{x) + Bq{x) , (72) 

where G(x) and q(x) are the gluon and quark densities, respectively, and A 
and B are the hard scattering cross sections calculable in perturbative QCD 
(pQCD). Similarly, the polarized cross section can be written as 

Aadi-jet(g, xP) = <J++.^^ - a+r^ = aAG(x) + 6Ag(x) , (73) 

where the first and second ± refer to the helicitics of the photon and nucleon, 
respectively. The double spin asymmetry for di-jet production is then 

^ ^ Ao-di-jct ^ a AG'(x) CTdMct ^ b Aq{x) 1 / o'd^ct j 
'^'^'"^ 2crdi-jot A G{x) 2(Tdi-jct B q(x) 2 i fJdi-jct / 
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The experimental asymmetry ^oxp in DIS is related to the photon asymmetry 

^cxp - PePNDAf-'^' , (75) 

where Pe and Pn are the electron and nucleon polarizations, respectively, and 
D is the depolarization factor of the photon. 

At low X, the gluon density dominates over the quark density, and thus the 
photon-gluon fusion subprocess dominates. There we simply have 

di-jct _ a AGjx) 

which provides a direct measurement of the gluon polarization. Because of the 
helicity selection rule, the photon and gluon must have opposite helicities to 
produce a quark and antiquark pair and hence a/ A = —I. Therefore, if AG{x) 
is positive, the spin asymmetry must be negative. Leading-order calculations 
in Refs. | 118| , 141, 245, HE] show that the asymmetry is large and is strongly 



sensitivitive to the gluon polarization. 

At NLO, the one-loop corrections for the PGF and QCDC subprocesses must 
be taken into account. In addition, three-jet events with two of the jets too close 
to be resolved must be treated as two-jet production. The sum of the virtual 
(2^2 processes with one loop) and real (2 — > 3 leading-order processes) 
corrections are independent of the infrared divergence. However, the two-jet 
cross section now depends on t he scheme in w hich the jets are defined. NLO 
calculations carried out in Refs. |233| , ^34 ^46| , show that the strong sensitivity 
of the cross section to the polarized gluon distribution survives. In terms of the 
spin asymmetry, the NLO effects do not significantly change the result. 

Since the invariant mass of the di-jet is itself a large mass scale, two-jet 
production can also be used to measure AG(x) even when the virtuality of 
the photon is small or zero (real photon). A great advantage of using nearly- 
real photons is that the cross section is large due to the infrared enhancement, 
and hence the statistics are high. An important disadvantage, however, is that 
there is now a contribution from the resolved photons. Because the photon 
is nearly on-shell, it has a complicated hadronic structure of its own. The 
structure can be described by quark and gluon distributions which have not 
yet been well determined experimentally. Some models of the spin-dependent 
parton distributions in the photon are discussed in Ref. [|155|. Leading-order 



calculations |27C, 9^ show that there are kinematic regions in which the resolved 
photon contribution is small and the experimental di-jet asymmetry can be used 
favorably to constrain the polarized gluon distribution. 



5.3 AG{x) from Large-p^ Hadron Production in e — p Scat- 
tering 

For e — p scattering at moderate center-of-mass energies, such as in fixed target 
experiments, jets are hard to identify because of their large angular spread and 
the low hadron multiplicity. However one still expects that the leading hadrons 
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in the final state reflect to a certain degree the original parton directions and 
flavors (discounting of course the transverse momentum, of order Aqcdj from 
the parton intrinsic motion in hadrons and from their fragmentation). If so, one 
can try to use the leading high-p^ hadrons to tag the partons produced in the 
hard subprocesses considered in the previous subsection. 

Bravar et al. |Q have proposed to use high-py hadrons to gain access to 
AG{x). To enhance the relative contribution from the photon-gluon fusion 
subprocess and hence the sensitivity of physical observables to the gluon distri- 
bution they propose a number of selection criteria for analysis of the data and 
then test these "cut" criteria in a Monte Carlo simulation of the COMPASS ex- 
periment. These simulations show that these cuts are effective in selecting the 
gluon-induced subprocess. Moreover, the spin asymmetry for the two-hadron 
production is large (10-20%) and is strongly sensitive to the gluon polarization. 

Because of the large invariant mass of the hadron pairs, the underlying 
subprocesses can still be described in perturbative QCD even if the virtuality of 



the photon is small or zero |144|. This enhances the data sample but introduces 
additional sub-processes to the high-p^ hadron production. The contribution 
from resolved photons, eg. from j ^ qq fluctuations, appears not to overwhelm 
the PGF contribution. Photons can also fluctuate into p mesons with p-nucleon 
scattering yielding large-p^ hadron pairs. Experimental information on this 
process can be used to subtract its contribution. After taking into account 
these contributions, it appears that the low-virtuality photons can be used as 
an effective probe of the gluon distribution to complement the data from DIS 
lepton scattering. 

5.4 AG{x) from Open-charm (Heavy-quark) Production in 
e — p Scattering 

Heavy quarks can be produced in e ~ p scattering through photon-gluon fusion 
and can be calculated in pQCD (see Fig. |l^). In the deep- inelastic scatter- 
ing region, the charm quark contribution to the gi{x, Q^) structure function is 



known |154| 



<?J(x,Q2) = ^ f'^Ap(^,QAAG{y,p'), (77) 

Jax V \y J 

where a = 1 + 4to^/Q^, and 

AF(x, ^ {2x - 1) In + 77(3 - 4a;) , (78) 

with 77^ = 1 — 4?7i^x/Q^(l — x). This result assumes that, because of the large 
charm quark mass, the direct charm contribution (eg. through Ac(x)) is small 
and the light-quark fragmentation production of charm mesons is suppressed. 
The X dependence of the structure function, if measured, can be deconvoluted 
to give the polarized gluon distribution. The renormalization scale p, can be 
taken to be twice the charm quark mass 2r7ic. 
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a) 



b) 



Figure 14: Feynman diagrams for charm production via Photon Gluon Fusion. 



Following Ref. [111|, the open charm electro-production cross section is large 
when is small or vanishes and can be written 



dQ'^dv 



where the virtual photon flux is 



ae,n 2(1 -y) + y^ + Qy2E^ 



2tt 



(79) 



(80) 



E and i' are the lepton and photon energies and y = v/E. For a fixed y, the flux 
is inversely proportional to Q^. The second factor in Eq. ^is the photonuclcon 
cross section. 

The cross section asymmetry is the simplest at the real-photon point = 0. 
The total parton cross section for photon-gluon fusion is 



87racmas(g) 
9i 



-/3(2-/32)-Hi(3-/34)lni±| 



(81) 



where (3 = -^/l — Am^/s is the center-of-mass velocity of the charm quark, and 
s = {q + xgP)^ is the invariant mass of the photon-gluon system. On the other 
hand, the spin-dependent cross section is 



Act 



9s 



3/3 - In 



l + f3' 
1-/3 



(82) 



The photon-nucleon asymmetry for open charm production can be obtained by 
convoluting the above cross sections with the gluon distribution, giving 



2MEy 
Aval 



dsAa{s)AGixG,s) 



r2MEy 
J Ami 



dsd{s)G(xG, s) 



(83) 
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Figure 15: Feynman diagrams for Direct photon production. 



where xq = s/2MpfEy is the gluon momentum fraction. Ignoring the depen- 
dence, the l—P spin asymmetry is related to the photon- nucleon spin asymmetry 
by A'[^ — DA'^'^j^, where D is the depolarization factor introduced before. 

The NLO corrections have recently been calculated by Bojak and Strat- 
mann |Q and Contogouris et al. |112|. The scale uncertainty is considerably 
reduced in NLO, but the dependence on the precise value of the charm quark 
mass is sizable at fixed target energies. 

Besides the total charm cross section, one can study the distributions of the 
cross section in the transverse momentum or rapidity of the charm quark. The 
benefit of doing this is that one can avoid the region of small xq where the 
asymmetry is very small ^7(| . 

Open charm production can be measured experimentally by detecting 13° 
mesons from charm quark fragmentation. On average, a charm quark has about 
60% probability of fragmenting into a D^. The D'^ meson can be reconstructed 
through its two-body decay mode 13° K~ + 7r+ ; the branching ratio is about 
4%. Additional background reduction can be achieved by tagging Z3°7r+ 
through detection of the additional 7r+ . 

J/iJj production is, in principle, also sensitive to the gluon densities. How- 
ever, because of ambiguities in the production mechanisms [184 , any informa- 
tion on AG{x) is likely to be highly model-dependent. 



5.5 AG{x) from Direct Photon Production in p — p Colli- 
sions 

AG{x) can be measured through direct (prompt) photon production in proton- 
proton or proton-antiproton scattering [Q. At tree level, the direct photon 
can be produced through two underlying subprocesses: Compton scattering 
qg — !■ 57 and quark-antiquark annihilation qq — s- 75, as shown in Fig. 
In proton-proton scattering, because the antiquark distribution is small, direct 
photon production is dominated by the Compton process and hence can be used 
to extract the gluon distribution directly. 
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Consider the collision of hadron A and B with momenta Pa and Pb , respec- 
tively. The invariant mass of the initial state is s = (Pa + Pb)^- Assume parton 
a (6) from the hadron A (B) carries longitudinal momentum XaPA [xaPs)- The 
Mandelstam variables for the parton subprocess a + 6 ^ 7 + c are 

S = XaXbS, t = Xat, U — XjjU , (84) 

where we have neglected the hadron mass. The parton-model cross section for 
inclusive direct-photon production is then 

^7$^ = E/ dxadx,f%{xa,^i')fB{x,,^^^)E,^ . (85) 

^ ab ^ 

For the polarized cross section Actab! the parton distributions fA.B are replaced 
by polarized distributions A/a^Bi and the parton cross sections aab are replaced 
by the spin-dependent cross section A&ab- The tree- level parton scattering cross 
section is 

E^-^ = ao^as^\M\^S{s + i+u) , (86) 

where the (5-function reduces the parton momentum integration into one inte- 
gration over, say, Xa with range [—u/{s + t),l] and 



For the polarized case, we have the same expression as in Eq. ( |8^ ) but with 



si ' ' ""^^^ 9 ici 

In the energy region where the Compton subprocess is dominant, we can 
write the proton-proton cross section in terms of the dccp-inelastic structure 



functions F2 and gi and the gluon distributions G and AG 70 



Ej '^^^^ = / dxadxb 
^'y~^^ = / dxadxb 



^^^^^^G{xb, fi')E,§^ + {Xa - Xb 

^J^l^^AG{x,,,W-^ + ^^a^x,) 



89) 



Here the factorization scale fi is usually taken as the photon transverse momen- 
tum pt- 

Unfortunately, the above simple picture of direct photon production is com- 
plicated by high-order QCD corrections. Starting at next-to- leading order the 
inclusive direct-photon production cross section is no longer well defined because 
of the infrared divergence arising when the photon momentum is coUincar with 
one of the final state partons. To absorb this divergence, an additional term 
must be added to Eq. ( ^5|) which represents the production of jets and their 
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subsequent fragmentation into photons. Therefore, the total photon production 
cross section depends also on these unknown parton-to-photon fragmentation 
functions. Moreover, the separation into direct photon and jet-fragmented pho- 
ton is scheme-dependent as the parton cross section E^ddab/d^P^ depends on 



the methods of infrared subtraction |147| 



To minimize the influence of the fragmentation contribution, one can impose 



an isolation cut on the experimental data |139|. Of course the parton cross 
section entering Eq. (|85| ) must be calculated in accordance with the cut criteria. 
An isolation cut has the additional benefit of excluding photons from vr*' or 77 
decay. When a high-energy 7r° decays, occasionally the two photons cannot 
be resolved in a detector or one of the photons may escape detection. These 
backgrounds usually reside in the cone of a jet and are largely excluded when 
an isolation cut is imposed. 

The NLO parton cross sections in direct photon production have been calcu- 



lated for both polarized and unpolarized scattering |147| . Comparison between 
the experimental data and theory for the latter case is still controversial. While 
the collider data at large pr are described very well by the NLO QCD calcula- 
tion H] , the fixed-target data and collider data at low-p^ are under-predicted by 
theory. Phenomenologically, this problem can be solved by introducing a broad- 
ening of the parton transverse momentum in the initial state [ pT| . Theoretical 
ideas attempting to resolve the discrepancy involve a resummation of threshold 



corrections [215| as well as a resummation of double logarithms involving the 



parton transverse momentum |212, 102|. Recently, it has been shown that a 



combination of both effects can reduce the discrepancy considerably [213| 



5.6 AG'(x) from Jet and Hadron Production 'm. p — p Col- 
lisions 

Jets are produced copiously in high-energy hadron colliders. The study of jets 
is now at a mature stage as the comparison between experimental data from 
Tevatron and other facilities and the NLO QCD calculations are in excellent 
agreement. Therefore, single and/or di-jet production in polarized colliders can 
be an excellent tool to measure the polarized parton distributions, particularly 
the gluon helicity distribution |Q . 

There are many underlying subprocesses contributing to leading-order jet 
production: qq' -> qq' , qq' -> qq' , qq qq, qq q'q', qq qq, qq gg, 
gg qq, qg qg, gg — > gg. Summing over all pairs of initial partons ab and 
subprocess channels ab cd, and folding in the parton distributions fa/Ai^a), 
etc., in the initial hadrons A and B, the net two-jet cross section is 

J dxadxbfa/A{xa)fb/B{xb)^p^{ab cd) . (90) 

abed 

For jets with large transverse momentum, it is clear that the valence quarks 
dominate the production. However, for intermediate and small transverse mo- 
mentum, jet production is overwhelmed by gluon-initiated subprocesses. 
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Studies of the NLO corrections are important in jet production because the 
QCD structure of the jets starts at this order. For polarized scattering, this has 



been investigated in a Monte Carlo simulation recently |113|. The main result of 
the study shows that the scale dependence is greatly reduced. Even though the 
jet asymmetry is small, because of the large abundance of jets, the statistical 
error is actually very small. 

Besides jets, one can also look for leading hadron production, just as in 
electroproduction considered previously. This is useful particularly when jet 
construction is difficult due to the limited geometrical coverage of the detectors. 
One generally expects that the hadron-production asymmetry has the same level 
of sensitivity to the gluon density as the jet asymmetry. 

5.7 Experimental Measurements 

The first information on AG has come from NLO fits to inclusive deep-inelastic 



scattering data as discussed in Sec. 3.5. Also recent semi-inclusive data from 



the HERMES experiment indicates a positive gluon polarization at a moderate 
XG- Future measurements from COMPASS at CERN, polarized RHIC, and 
polarized HERA promise to provide much more accurate data. 

5.7.1 Inclusive DIS Scattering 



As discussed in Sec. 3.5, the spin-dependent structure function Q^) is sen- 
sitive to the gluon distribution at NLO. However, to extract the gluon dis- 
tribution, which appears as an additive term, one relies on the different Q^- 
dependence of the quark and gluon contributions. 

The biggest uncertainty in the procedure of the NLO fits is the parametric 
form of the gluon distribution at Qq- It is known that by taking different 
parameterizations, one can get quite different results. 

5.7.2 HERMES Semi-inclusive Scattering 



The HERMES experiment has been described in Sec. 2.3, In a recent pub- 
lication [Q, the HERMES collaboration reported a first measurement of the 
longitudinal spin asymmetry = —0.28 ± 0.12 ± 0.02 in the photoproduction 
of pairs of hadrons with high transverse momentum p^, which translate into a 
(AG/G) = 0.41 ± 0.18 ± 0.03 at an average {xg) = 0.17. 

Following the proposal of Ref. js^, the data sample contains hadron pairs 
with opposite electric charge. The momentum of the hadron is required to be 
above 4.5 GeV/c with a transverse component above 0.5 GeV/c. The minimum 
value of the invariant mass of the two hadrons, in the case of two pions, is 1.0 
GeV/c^. A nonzero asymmetry is observed if the pairs with > 1.5 GeV/c 
andp^" > l.OGeV/c are selected. The measured asymmetry is shown in Fig. ^ 
with an average of 0.06(GeV/c)^. If p'^^ > 1.5 GeV/c is not enforced the 
asymmetry is consistent with zero. 
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Figure 16: Data for high pT hadrons from the HERMES experiment 
(left) the asymmetry and (right) the extracted gluon polarization. 



52| showing 



The measured asymmetry was interpreted in terms of the following processes: 
lowest-order deep-inelastic scattering, vector-dominance of the photon, resolved 
photon, and hard QCD processes - Photon Gluon Fusion and QCD Compton 



effects. The PYTHIA [263 Monte Carlo generator was used to provide a model 
for the data. In the region of phase space where a negative asymmetry is 
observed, the simulated cross section is dominated by photon gluon fusion. The 
sensitivity of the measured asymmetry to the polarized gluon distribution is 
also shown Fig. Note that the analysis does not include NLO contributions 
which could be important. The HERMES collaboration will have more data on 
this process in the near future. 



5.7.3 COMPASS Experiment 

The COMPASS expriment at CERN will use a high-energy (up to 200 GeV) 
muon beam to perform deep-inelastic scattering on nucleon targets, detecting 
final state hadron production The main goal of the experiment is to 



measure the cross section asymmetry for open charm production to extract the 
gluon polarization AG. 

For the charm production process, COMPASS estimates a charm production 
cross section of 200 to 350 nb. With a luminosity of 4.3 x 10'^^ cm~^day"^, they 
predict about 82,000 charm events in this kinematic region per day. Taking into 
account branching ratios, the geometrical acceptance and target rescattering, 
etc., 900 of these events can be reconstructed per day. The number of back- 
ground events is on the order of 3000 per day. Therefore the total statistical 
error on the spin asymmetry will be about 5A'i^j^ — 0.076. 

Shown in Fig. |l^ are the predicted asymmetries A'ji^j^ and A'^^^ for open 
charm production as a function of y. The curves correspond to three differ- 
ent models for AG. From the results at different y, one hopes to get some 
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Figure 17: Predictions of the open charm asymmetry for the COMPASS exper- 
iment. The curves are predictions for three representative gluon spin distribu- 
tions 1 151 1 . Also shown are typical error bars expected from the measurement. 



information about the variation of AG as a function of x. Measurements with 
high -pT hadrons will also be used to complement the information from charm 
production. 



5.7.4 l\G{x) from RHIC Spin Experiments 

One of the primary goals of the RHIC spin experiments is to determine the 
polarized gluon distribution. This can be done with direct photon, jet, and 
heavy quark production. 

Direct photon production is unique at RHIC. This can either be done on 
inclusive direct photon events (PHENIX) or photon-plus-jet events (STAR). 
Estimates of the background from qq annihilation show a small effect. Shown in 
Fig. |l| is the sensitivity of STAR measurements of AG(a;) in the channel 
7-|-jet-|-X. The solid line is the input distribution and the data points represent 
the reconstructed AG'(x). For inclusive direct photon events, simulations show 
very different spin asymmetries from different spin-dependent gluon densities. 

Jet and heavy flavor productions are also favorable channels to measure 
polarized gluons at RHIC. The interested reader can consult the recent review 
in Ref. 



5.7.5 AG(x) from Polarized HERA 

The idea of a polarized HERA collider (e~ — p) has been described in Sec 2.4. 
Here we highlight a few experiments which can provide a good measurement of 
the polarized gluon distribution |12C|. 

First of all, polarized HERA will provide access to very large and low 
X regions compared with fixed-target experiments. At large x, can be as 
large as lO** GeV^. Thus, one can probe the gluon distribution through the 
variation of the gi structure function. An estimate from an NLO pQCD analysis 
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Figure 18: Projected sensitivity 147 to the gluon polarization from direct pho- 
ton production at RHIC with the PHENIX detector. The curves represent 



different assumptions for the gluon spin distribution 151 156 



shows that the polarized HERA data on gi can reduce the uncertainty on the 
total gluon helicity to ±0.2(exp)±0.3(theory). 

Polarized HERA can also measure the polarized gluon distribution through 
di-jet production. Assuming a luminosity 500 pb~^ and with the event selection 
criteria 5 < < 100 GeV^, 0.3 < y < 0.85 and > 5 GeV, the expected 
error bars on the extracted AG{x) are shown in Fig. The measured x region 
covers 0.002 < xq < 0.2. AG{x) can also be measured at polarized HERA 
through high-pT hadrons and jet production with real photons. 



6 Transverse Spin Physics 

6.1 The g2(x,Q^) Structure Function of the Nucleon 



As discussed in Sec. |3.2| , the structure function g2{x,Q'^) can be measured 
with a longitudinally polarized lepton beam incident on a transversely polar- 
ized nucleon target. For many years, theorists have searched for a physical 
interpr etation o f g2{x,Q^) in terms of a generalization of Feynman's parton 



model 142, 210| , as most of the known high-energy processes can be understood 



in terms of mco/ierenf scattering of massless, on-shell and collinear partons | 142 |. 
It turns out, however, that g2 is an example of a higher-twist structure function. 
Higher-twist processes cannot be understood in terms of the simple parton 



model [138|. Instead, one has to consider parton correlations initially present 
in the participating hadrons. Higher- twist processes can be described in terms 
of coherent parton scattering in the sense that more than one parton from a 
particular hadron takes part in the scattering. Higher-twist observables are in- 
teresting because they represent the quark and gluon correlations in the nucleon 
which cannot otherwise be studied. 



48 



Figure 19: Projected statistical and systematic error bars for the polarized 
gluon distribution using two-jet production with polarized e — p collisions at 
HERA. The curves represent different assumptions about the gluon spin distri- 
bution ]l51|. 



Why does g2{x, Q^) contain information about quark and gluon correlations? 
According to the optical theorem, 92(2;, Q^) is the imaginary part of the spin- 
dependent Compton amplitude for the process 

7*(+l) + A^(l/2) ^7*(0) + A^(-l/2) (91) 

where 7* and N represent the virtual photon and nucleon, respectively, and 
the labels in the brackets are helicities. Thus Compton scattering involves a 
t-channel helicity exchange. When the process is factorized in terms of par- 
ton subprocesses, the intermediate partons must carry this helicity exchange. 
Because of the vector coupling, massless quarks cannot undergo helicity flip in 
perturbative processes. Nonetheless, the required helicity exchange is fulfilled 
in two ways in QCD: first, through single quark scattering in which the quark 
carries one unit of orbital angular momentum through its transverse momen- 
tum; second, through quark scattering with an additional transversely-polarized 
gluon from the nucleon target. These two mechanisms are combined in such a 
way to yield a gauge- invariant result. 

To leading order in a^, g2{x B,Q^) can be expressed in terms of a simple 
parton distribution Agy (a;) [|l72| , 

52(0:5, Q2) = i ^ e2(Ag,T(a;B, Q') + /^q^T{xB, Q')) , (92) 

■i 

where 

A<ZT(:r) = ^/ ^e'^^-(P5x|g-(0)7x75<z(An)|P5^) , (93) 
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and S± is the transverse polarization vector and j± is the component along the 



same direction. Although it allows a simple estimate of 32 in the models |180|, 



the above expression is deceptive in its physical content. It has led to incorrect 



identifications of twist-three operators pq, 258 1 and incorrect next-to-leading 



order coefficient functions 1 204 1 . When the leading- logarithmic corrections were 
studied, it was found that AqT{x, Q^) mixes with other distributions under scale 
evolution |261|. In fact, Agr(x,Q^) is a special moment of more general parton 



distributions involving two light-cone variables 



Aqrix)^- f dy{Ki{x,y)+K2{x,y)) , (94) 



X 



where the Ki{x,y) are defined as 

^^^^x\+^^.{y-x) (^ps\Tl,(o)iD" {fin)ij{Xn)\PS) 
= S''-f5^Ki{x,y)+iT'' ^K2ix,y) + ... (95) 
where T" = e"^''^ S±Pjns. Under a scale transformation, the general distribu- 



tions Ki{x,y) evolve autonomously while the AqiT{x) do not l9d, 131]. The 



first result for the leading logarithmic evolution of the twist-three distributions 



(and operators) [Q has now been confirmed by many studies |252, Q. 

Thus an all-order 52 factorization formula is much more subtle than is in- 
dicated by the leading-order result. It involves the generalized two- variable 
distributions, Ki{x,y), 

9t{xb,Q'^) [ ^^^^^ (Ca( — ,—,as] Ki{x,y) + {xB -xb)] , 

'^J-i xy \ \x y J J 

(96) 

where Ca are the coefficient functions with a summing over different quark 
flavors and over gluons. Accordingly, a perturbative calculation of g2 in terms 



of quark and gluon external states must be interpreted carefully 1 22S | . Recently, 



the complete one-loop radiative corrections to the singlet and non-singlet 32 have 



been published |192|. The result is represented as the order-ofg term in Ci and 
is one of the necessary ingredients for a NLO analysis of 52 data. Note that 
the Burkhardt-Cottingham sum rule, 32(2^, Q^)dx — 0, survives the radiative 
corrections provided the order of integrations can be exchanged . 

As an example of the interesting physics associated with 172, we consider its 



second moment in 



X 



[' dxxg2{x,Q^)^l{~a2{Q^) + d2{Q^)) , (97) 
Jo •J 

where a2{Q^) is the second moment of the gi{x) structure function. Here d2{Q^) 
is the matrix element of a twist-three operator, 

{PS\^i)gF''^^'-^''^\PS) = 2d2S^'' P^^'^ P"^ , (98) 
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where = {l/2)e'"'°''^ Fa/i, and the different brackets -(•••) and [•••]- denote 
symmetrization and antisymmetrization of indices, respectively. The structure 
of this twist-three operator suggests that it measures a quark and a gluon am- 
pHtude in the initial nucleon wave function. 

To better understand the significance of d2{Q^), we consider a polarized 
nucleon in its rest frame and consider how the gluon field inside of the nucleon 
responds to the polarization. Intuitively, because of parity conservation, the 
color magnetic field B can be induced along the nucleon polarization and the 
color electric field E in the plane perpendicular to the polarization. Introducing 
the color-singlet operators Ob — ip'^gBip and Oe = 'ip^d x gEip, we define the 
gluon-field polarizabilities xb and xe in the rest frame of the nucleon, 

(PSIObMPS) = xb,e2M^S . (99) 

Then it is easy to show 

^2 = (2xb+Xb)/3 . (100) 

Thus ^2 measures the response of the color electric and magnetic fields to the 
polarization of the nucleon. 

The experimental measurements of the 52 structure function started with 
U) and E142jp[ ||l collaborations. Subsequently, the E143 g, 



the SMC 



E154 I?!, and E155 collaborations have also measured and published 

their data. The combined E143 and E155 data for proton and deutero n ar e 
shown in Fig. The solid line shows the twist- two contribution to 52 only |278l. 



The dashed and dash-dotted lines are the bag model calculations by Song [ 267 | 



and Stratmann |269 



Neglecting the contributions from x < 0.02 and x > 0.8 and the de- 
pendence, the E155 collaboration Q has integrated their data to get 
/ dxg^ = -0.022 ± 0.071 and / dxgf= 0.023 ± 0.044. The resuhs are con- 
sistent with the Burkardt-Cottingham sum rules within the relatively large 
errors. The second moments allow an extraction of the d2 matrix elements. 
E155 found = 0.005 ± 0.008 and d'i = 0.008 ± 0.005 at an average of 
5 GeV^. A combined analysis of the E142, E143, E154, and E155 data yields 
dP = 0.007± 0.004 and d^ = 0.004± 0.010. These numbers are generally consis- 



tent with bag model |267, 269, 19S , |193[ and chiral quark model [279| estimates, 
and are 1 to 2a away from QCD sum rule calculations | p68| , |6^ , 134 1. The error 
bars on the present lattice calculation are still relatively large [ |159[ . 

According to the simple quark model, the d2 matrix element in the neutron 
should be much smaller than that in the proton because of SU(6) spin-fiavor 
symmetry. While the proton ^2 has been constrained with reasonable precision, 
the neutron c?2 has a much larger error bar. In the near future, JLab experiments 
with a polarized '^He target [59| can improve the present error on the neutron 
^2 and hence test the quark model predictions. 

6.2 Tranversity Distribution 

Along with the unpolarized and polarized quark distributions - qi{x,Q'^) and 
Aqi{x, Q^) - discussed above, a third quark distribution exists at the same order 
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Figure 20: (Upper) Data for 52 for the proton and deuteron from E143 |^ and 
E155 collaborations. (Lower) Comparison between data and calculations 
for d2 , the second moment of 32 ■ 
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(twist two) as the other two distributions. Note that no corresponding transverse 
spin distribution exists for gluons (due to hehcity conservation). 

This transversity distribution, 5qi(x,Q'^), can be described in the Quark- 
Parton Model as the difference in the distribution of quarks with spin aligned 
along the nucleon spin vs. anti-aligned for a nucleon polarized transverse to its 
momentum. The structure function related to transversity is given by 



/ii(x,g2) 



e-%(a;,Q^) 



(101) 



The first moment of the transversity distributions also leads to an interesting 
observable - the nuclcon's tensor charge a\: 



[5q^{x) - 5qi{x)] dx . 



(102) 



In terms of nucleon matrix elements, this tensor charge is defined |16S] as: 

2a\ {S'^P" - S^P") = {PS\'4j,a'"'ij5^^\PS) . (103) 

Recent calculations have made estimates of the tensor charges using QCD Sum 
Rules [|l69[ |l9||, Lattice QCD and within the Chiral Quark Model ]202 |. 

In a non-relativistic model the transversity is equal to the longitudinal spin 
distribution {Sqi = Aq.;) because the distribution would be invariant under 
the combination of a rotation and a Lorentz boost. Relativistically, this is 
not the case and 5qi could be significantly different from Aqi. The challenge 
to gaining experimental information on Sqi lies in its chiral structure. In the 
helicity basis |181, 182 1 Sqi represents a quark helicity flip, which cannot occur 
in any hard process for massless quarks within QED or QCD. This chiral-odd 
property of transversity makes it unobservable in inclusive DIS. In order to 
observe hi{x, Q^) a second non-perturbative process that is a lso c hiral-odd must 
take place. This was first discussed by Ralston and Soffer |251| in connection 
witli Drell-Yan production of di-muons in polarized p — p collisions. Here the 
transversity distribution of both protons results in a chiral-even interaction. 

Several calculations have suggested that the transversity distribution may 



be accessible in semi-inclusive lepton-nucleon scattering |108, 183, 20(;, 236, 
In this process a chiral-odd fragmentation function, leading to a lepto- 
produced hadron, offsets the chiral-odd transversity distribution. Many of these 
calculations take advantage of an inequality 



\Sq,ix)\ < 



q.i{x) + Aqi{x) 



(104) 



discovered by Soffer ]264[| to limit the possible magnitude of hi{x). 



Calculations |206, 23(:, 207] have also detailed a set of spin distribution and 
fragmentation functions that are accessible from leading and next-to-leading 
twist processes. In fact in some cases the next-to-leading twist processes can 
dominate, especially at low and with longitudinally polarized targets. 
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Figure 21: Single spin asymmetry for pion electroproduction from the HERMES 
experiment 133] vs. pr (left) and x (right). 



Potentially relevant experimental information has recently come from the 
HERMES collaboration and SMC collaboration. HERMES has measured the 
single-spin azimuthal asymmetry for pions produced in deep-inelastic scattering 
of unpolarized positrons from a longitudinally polarized hydrogen target p3[ 
(see Fig. ^l|). A related measurement has been reported by SMC using 
a transversely polarized target. The HERMES asymmetry is consistent with 
a sint/) distribution, where is the angle between the lepton scattering plane 
and the plane formed by the virtual photon and pion momenta. While the x 
dependence of the asymmetry is relatively weak except for the smallest x = 0.04 
point, the pr dependence shows a rapid rise up to 0.8 GeV~^. The average 7r+ 
asymmetry averaged over the full acceptance is 0.022 ± 0.005 ± 0.003 while 
the asymmetry for 7r~ production is consistent with zero. Some models [ 133 | 
suggest that the product of the transversity distribution times the chiral-odd 
fragmentation function can account for the observed asymmetry. 

Transversity may also play a role in observed single-spin asymmetries inp—p 
collisions. These possibilities are discussed in the next section. 



6.3 Single-Spin Asymmetries Prom Strong Interactions 

As we have discussed in Sec. |4.2| , single-spin asymmetries can arise from pro- 
cesses involving parity- violating interactions, such as W and Z boson production 
in p — p collisions at RHIC. In this subsection, we discuss a different class of 
single-spin asymmetries which are generated entirely from strong interaction ef- 
fects. While we do not have enough space here to make a thorough examination 
of the subject, we briefly discuss the phenomena, a few leading theoretical ideas, 
and some other related topics. A recent review of the subject can be found in 
Ref. I222I. 
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For more than two decades, it has been known that in hadron-hadron scat- 
tering with one beam transversely polarized, the single-particle inclusive yield 
at nonzero pt has an azimuthal dependence in a coordinate system where z is 
chosen to lie along the direction of the polarized beam, and x along the beam 
polarization | |203{ . It is easy to see that the angular dependence is allowed by 
strong interaction dynamics. If the momentum of the polarized beam is pt, and 
that of the observed particle p'o, the angular distribution reflects the existence 
of a triple correlation, 

Pb^Po- S , (105) 

where S is the beam polarization. The correlation conserves parity and hence 
is not forbidden in strong interactions. Although it is nominally time-reversal 
odd, the minus sign can be canceled, under the time-reversal transformation, by 
a factor of i from an interference of two amplitudes with different phase factors. 
The angular correlation is usually characterized by the spin asymmetry 

where da"^'^ are the cross sections with reversed polarizations, and pj- is the 
transverse momentum of the produced particle, xp is the Feynman x variable, 
xp — Pl/p'l'^^j where is the longitudinal momentum of the produced hadron 
and p'^'^^ is the maximum allowed longitudinal momentum. An example of a 
single spin asymmetry for tt production is shown in Fig. p2. After examining 



the existing data, one finds the following interesting systematic effects [222\: 



A^v is significant only in the fragmentation region of the polarized beam. 
It increases almost linearly with xp when the target is unpolarized. 

Ajv is large only for moderate transverse momentum px ■ 

Am and its sign show a strong dependence on the type of polarized beam 
{p,p) and produced particles (tt^jTTo). 



That is strikingly large is the most impressive aspect of the phenomenon. 

The simplest theory explaining Ajsr is one that assumes an underlying parton 
process: partons from the parent hadrons scatter and fragment to produce the 
observed particle. To get the single spin asymmetry, one requires, for instance, 
that quarks change their helicity during hard scattering. However, chiral sym- 
metry then dictates that the asymmetry is proportional to the quark mass mq 
which is vanishingly small for light quarks. Thus the simple parton model for 



An cannot yield the magnitude of the observed symmetry |20C|. 

For the moment, the leading theoretical ideas in the literature are still based 
on the parton degrees of freedom. However, the spin-flip is introduced through 
more complicated mechanisms: Either the initial and final partons are assumed 
to have novel nonpertrubative distribution and fragmentation functions, respec- 
tively, or the parton hard scattering involves coherent processes. 

In the latter case, the asymmetry can arise from the coupling of chiral even 
(odd) twist-two (twist-three) parton correlations in the polarized nuclcon and 
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Figure 22: Measured single-spin asymmetries from tt^ production in p — p scat- 
tering (left) and i:^ production in p ~ p and p~p scattering (right) at Fermi- 
lab ||203|. The curves are theoretical fits from Ref. [244|. 



chiral even (o dd) twist-three (tw ist-two) fragmentation functions of the scat- 



tered partons ]132| , 243, 185, 244 1. The required phase difference is generated 
from the interference of the hard scattering amplitudes in which one of the hard 
propagators is on-shell. The predicted asymmetry is of order Aqcd /pt in the 
large pt limit, which is a characteristic twist-three effect. For moderate px, 
Ajsf can be a slowly decreasing function of px |244|. The comparison between 
the experimental data and the phenomenological prediction seems to yield good 
agreement [244|. It is not clear, however, that the available fixed target data can 
be fully described by perturbative parton scattering. One needs more data at 
higher energy to test the scaling property inherent in a perturbative description. 

The alternative is to consider nonperturbative mechanisms to generate the 
phase difference. This can be done by introducing transverse-momentum depen- 
dent parton distributions |262| and fragmentation functions [108 . In a trans- 



versely polarized nucleon, the transverse momentum distribution may not be 
rotationally invariant. It may depend on the relative orientation of the spin and 
momentum vectors. Likewise, when a transversely polarized quark fragments, 
the amplitude for hadron production can depend on the relative orientation be- 
tween the hadron momentum and the quark spin. Both mechanisms have been 
shown to produce large single spin asymmetries @, [t^, |7^, ^ . Here again the 
applicability of the model for the existing data is not clear. In particular, the 
fitted fragmentation functions and parton distributions must be tested in dif- 
ferent kinematic regions. Moreover, the new distributions do not possess color 
gauge invariance. 

A phenomenological model for parton scattering with formation of large-pT 
hadrons was proposed by Boros, Liang and Meng Although not derived 

from field theory, the model has a very intuitive physical picture and successfully 
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describes the data. It would be interesting to test the predictive power of the 
model in future experiments. The RHIC spin facility can test many of these 
theoretical ideas with a variety of experimental probes including |254| polarized 
Drell-Yan, dimeson production, etc. 

A subject closely related to the single-spin asymmetry is the polarization 



of hyperons, such as A, produced in unpolarized hadron collisions [171|. The 
observed polarization is perpendicular to the plane formed by the beam and 
hyperon momenta. Many theoretical models have been invented to explain the 
polarization |l40| . Most models are closely related to those devised to explain 
the single spin asymmetry. 



7 Off-Forward Parton Distributions 



In this section, we discuss some of the recent theoretical developments on gen- 
eralized (off-forward) parton distributions (OFPD) and their relation to the 
angular momentum distributions in the nucleon. We will also consider possible 
experimental processes, such as deeply virtual Compton scattering (DVCS) and 
meson production, to measure these novel distributions. 

OFPD's were first introduced in Ref. |123{ and discovered independently 
in Ref. |186| in studying the spin structure of the nucleon. Radyushkin and 
others have introduced slightl y di f ferent ver s ions of the distributions, but the 
physical content is the same [^47|, 248 24£, 11C|. The other names for these 



functions range from off-diagonal, non-forward and skewed to generalized parton 
distributions. Here we follow the discussion in Ref. |189|. 

One of the most important sources of information about the nucleon struc- 
ture is the form factors of the electroweak currents. It is well known that the 
vector current yields two form factors 



_ irrf^'^n 



(107) 



where = P' — P and Fi and F2 are the Dirac and Pauli form factors, respec- 
tively. F2 gives the anomalous magnetic moment of the nucleon, k = i^2(0). 
The charge radius of the nucleon is defined by 



= -6 



(108) 



with Ge = Fi- Q^/{U'P)F2 
factors, 



The axial vector current also defines two form 



(P'|^7'^75^/^|P) = G^(Q^)C/7^75C/ + Gp{Q'')U^UiP) 



(109) 



The axial form factor Ga is related to the fraction of the nucleon spin carried by 
the spin of the quarks, AS, and can be measured from polarized deep inelastic 
scattering as discussed in previous sections. The pseudoscalar charge, Gp{0), 
can be measured in muon capture. 
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A generalization of the electroweak currents can be made through the fol- 
lowing sets of twist-two operators, 

O^i-''- = V'q7^^'7.3«-D^'---»£'^"V, (110) 

where all indices Hi - ■■ are symmetric and traceless as indicated by (...) in the 
superscripts. These operators form the totally symmetric representation of the 
Lorentz group. One can also introduce gluon currents through the operators: 

O^i-i^n = F'^'^'"iD^---iD'^^-^FJ^"K (111) 

For n > 1, the above operators are not conserved currents from any global 
symmetry. Consequently, their matrix elements depend on the momentum- 
transfer scale n at which they are probed. For the same reason, there is no 
low-energy probe that couples to these currents. 

One can then define the generalized charges a„(/i^) from the forward matrix 
elements of these currents 

(P|OMi-A*2|p^ = 2an{n^)P^'''P''' ■■■Pf'-K (112) 

The moments of the Feynman parton distribution q{x,fi^) are related to these 

charges 

j ^ dxx-^-^qix, i/) = ^ dxx^-' [q{x, /x') + (-l)"5(a;, m')] = a„(M') , (113) 

where q{x, fx^) is defined in the range — 1 < x < 1. For x > 0, q{x, fx^) is just the 
density of quarks which carry the fraction x of the parent micleon momentum. 
The density of antiquarks is customarily denoted as q(x, jj?), which in the above 
notation is — g(— a;,/x^) for x < 0. 

One can also define the form factors (Aq„.„i(t), (f), and Cqn{t)) of 

these currents using constraints from charge conjugation, parity, time-reversal 
and Lorentz symmetries 

(P'|0^i-'^"|P) 

= U{P')^^<''U{P) Aqnaiit)^"^ ■ ■ ■ A>'''+^P^''+' ■ ■ ■ P^"' 

i=0 

+ U{P' f^''^2^- U{P) Y B,„,2i(i)A''^...A''-+^P''^*+^..p'"'' 

1=0 

+ Cgn{t)Mod{n + 1, 2) -^C7(P')f/(P)A('^^ • • • A"-) , (114) 

where U{P') and U{P) are Dirac spinors, A^ = (P' -Pf =t,P= (P' + P)/2 
and Mod(n -|- 1, 2) is 1 when n is even and when n is odd. Thus Cg„ is present 
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only when n is even. We suppress the renormahzation scale dependence for sim- 
plicity. In high energy experiments, it is difficult to isolate the individual form 
factors. Instead it is useful to consolidate them into generalized distributions — 
the off- forward parton distributions (OFPD's). To accomplish this a light-light 
vector n'^ (n^ = 0) is chosen such that 

n-P = l , C = -n- A/2. (115) 

Then, 

n^,-n^AP'\0^'-''"\P)=H„{U)U jlU + E^UW " £/ , (116) 

where i/„(^,i) and En{^,t) are polynomials in of degree n/2 (n even) or 
n — 1/2 {n odd). The coefficients of the polynomials are form factors. The 
OFPD E{x, ^, t) and H{x, ^, t) are then defined as: 

J dxx"-i£(x,C,i) = En{^,t) 

j ^dxx''-^H{x,^,t) = Hn{^,t). (117) 

Since all form factors are real, the new distributions are also real. Moreover, 
because of time-reversal and hermiticity, they are even functions of ^. 

The OFPD's are more complicated than the Feynman parton distributions 
because of their dependence on the momentum transfer A. As such, they contain 
two more scalar variables besides the x variable. The variable t is the usual t- 
channel invariant which is always present in a form factor. The ^ variable 
is a natural product of marrying the concepts of the parton distribution and 
the form factor: The former requires the presence of a prefered momentum 
along which the partons are predominantly moving, and the latter requires a 
four-momentum transfer A; ^ is just a scalar product of these two momenta. 

7.1 Properties of the OfF-ForwEird Pcirton Distributions 

The physical interpretation of parton distributions is transparent only in light- 
cone coordinates and light-cone gauge. To sec this, wc sum up all the local 
twist-two operators into a light-cone bilocal operator and express the parton 
distributions in terms of the latter, 

=H,{x,u) Win mP)+E,{x,U)mP')'^^^U{P). 

(118) 

The light-cone bilocal operator (or light-ray operator) arises frequently in hard 
scattering processes in which partons propagate along the light-cone. In the 



(-fn) rlPe 



f\/2 



da n-A{an) 



V, (in) 
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light-cone gauge n ■ A = 0, the gauge link between the quark fields can be 
ignored. Using the light-cone coordinate system 



1 



(119) 



we can expand the Dirac field 
, _ , , dk+d'^ki 

,XJ_) 



2fc+(27r)3 
d{{k+,k^)v),{k)e 



e{k+)J2 (&A(fc+,fc±)MA(fc)e-^(" k+-S^-k^) 
\=± 

{x~ —x± ■k± ) 



(120) 



where V'-i- ~ P+'>P f^nd P± = ^7^7 . The quark (antiquark) creation and 
annihilation operators, (d^j.) and b\k {d\k), obey the usual commutation 
relation. Substituting the above into Eq. ( |118| ), we have |18£| 

, If d^k\ 

Fgix,0^ 

P' 
P' 



2p+V J 2^|a;2-^2|(2^)3 ^ 
b\ f(x-f)p+,fc| +A 



X < 



l((x-Op+,fcj 
dx (^{-x + Op+,-k± 



A 

b\ ([X 



((x + C)p+,fc± 
AAb^J{x + Op 



P 



P' 



dx(i-x + ^)p+,k^ 



for a; > ^ 
for ^ > X > 
for X < — ^ 

(121) 



where 1^ is a volume factor. The distribution has different physical interpre- 
tations in the three different regions. In the region x > ^, it is the amplitude 
for taking a quark of momentum k out of the nucleon, changing its momen- 
tum to fc -f A, and inserting it back to form a recoiled nucleon. In the region 
^ > X > — ^, it is the amplitude for taking out a quark and antiquark pair with 
momentum —A. Finally, in the region x < — ^, we have the same situation as 
in the first, except the quark is replaced by an antiquark. The first and third 
regions are similar to those present in ordinary parton distributions, while the 
middle region is similar to that in a meson amplitude. 

By recalling the definition of Jq g(fi) in terms of the QCD energy-momentum 



tensor T^^g 



,1 



P- 



d^x{x X Tq^g), 



1 



p- 



(122) 



it is clear that they can be extracted from the form factors of the quark and 
gluon parts of the T^^. Specializing Eq. (114) to (n = 2), 

{P'\mP) - U{P')Wg{t)^^''P"^ +Bg,^g{t)P^''ia-^''A^/2M 



+ Cq^g{t)/\^''A'''^ /m]u{P). 



(123) 
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Taking the forward limit of the ^ — component and integrating over three- 
space, one finds that the A, ^(0) give the momentum fractions of the nucleon 
carried by quarks and gluons {Aq{0) + ^g(O) = 1). On the other hand, substi- 



tuting the above into the nucleon matrix element of Eq. (122), one finds 118 



J,,9^7:[A,A0) + Bg,gm- (124) 



Therefore, the matrix elements of the energy-momentum tensor provide the 
fractions of the nucleon spin carried by quarks and gluons. There is an analogy 
for this. If one knows the Dirac and Pauli form factors of the electromagnetic 
current, Fi{Q'^) and F2{Q^), the magnetic moment of the nucleon, defined as 
the matrix element of (1/2)/ ^^^(x x JY, is Fi(0) + ^2(0). 

Since the quark and gluon energy-momentum tensors are just the twist- 
two, spin-two, parton helicity-independent operators, we immediately have the 
following sum rule from the off-forward distributions; 

1 

dxx[Hg{x, t) + Eg{x, ^, t)] = Ag{t) + Bg{t) , (125) 

-1 

where the ^ dependence, or Cq(t) contamination, drops out. Extrapolating the 
sum rule to t = 0, the total quark contribution to the nucleon spin is obtained. 
When combined with measurements of the quark spin contribution via polarized 
DIS measurements, the quark orbital contribution to the nucleon spin can be 
extracted. A similar sum rule exists for gluons. Thus a deep understanding of 
the spin structure of the nucleon may be achieved by measuring OFPD's in high 
energy experiments. 

A few rigorous results about OFPD's are known. First of all, in the limit 
^ — > and i — > 0, they reduce to the ordinary parton distributions. For instance, 

Hq{x,0,0) = q{x), 

Hg{x,0,0) = Aq{x), (126) 

where q(x) and Aq{x) are the unpolarized and polarized quark densities. Similar 
equations hold for gluon distributions. For practical purposes, in the kinematic 
region where 

y]<|«A/Ar and ^ « a; (127) 

an off-forward distribution may be approximated by the corresponding forward 
one. The first condition, ^y\t\ <<: M^, is crucial — otherwise there is a significant 
form-factor suppression which cannot be neglected at any x and ^. For a given 
^ is restricted to 

1^1 < - t/A) . (128) 

Therefore, when yj\t\ is small, ^ is automatically limited and there is in fact a 
large region of x where the forward approximation holds. 

The first moments of the off-forward distributions are constrained by the 
form factors of the electromagnetic and axial currents. Indeed, by integrating 
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over a;, we have [186] 



Flit) ,J ^dxEg{x,^,t)=F^{t) , 
»i 



dxHg{x,^,t) = G\{t) J dxEq{x,S„t) = G%{t) 



(129) 



-1 



where Fi, F2, Ga and Gp are the Dirac, Pauh, axial, and pseudo-scalar elastic 
form factors, respectively. The t dependence of the form factors are character- 
ized by hadron mass scales. Therefore, it is reasonable to speculate that similar 
mass scales control the t dependence of the off-forward distributions. 



The first calculation of the OFPD has been done in the MIT bag model |194|. 
The parameters are adjusted so that the electromagnetic form factors and the 
Feynman parton distributions are well reproduced. The shapes of the distribu- 
tions as a function of x are rather similar at different t and f . The t dependence 
of the energy-momentum form factors is controlled by a mass parameter be- 
tween 0.5 and 1 GeV^. The same distributions were also studied in the chiral 



quark-soliton model by Petrov et al. |242|. In contrast to the bag model results, 
the chiral soliton model yields a rather strong ^ dependence. The model also 
predicts qualitatively different behaviors in the regions |a:| > ^ and < ^, 
in line with the physical interpretation of the distributions. In the case of the 



distribution, the pion pole contribution is important [14S|. The OFPD's 
have als o be en modeled directly without a theory of the structure of the nucleon. 
In Ref . 1 274 , the distributions are assumed to be a product of the usual parton 
distributions and some ^-dependent form factors, independent of the variable ^. 
In Ref. |24£, 25C|, the so-called double distributions are modeled in a similar 
ansatz from which a strong ^ dependence is generated. 

Scale evolution of the OFPD's has received wide attention and is now com- 
pletely solved up to two loops. In the operator form, the evolution has been 
studied at the leading logarithmic approximation long before U227| . In terms of 
the actual distributions, the evolution equations at the leading-log can be found 
in Refs. ||l2|, |l8^, |4^, |24|, |24|, |lO| in different cases and forms. 

In a series of interesting papers, Belitsky and Miiller have calculated the 
evolution of the off- forward distributions at two loops |69) . The key observation 
is that perturbative QCD is approximately conformally invariant. The breaking 
of the conformal symmetry can be studied through conformal Ward identities, 
which allows one to obtain the two-loop anomalous dimension by calculating 
just the one-loop conformal anomaly. 



7.2 Deeply Virtual Exclusive Scattering 

Of course the eventual utility of the OFPD's depends on whether they can ac- 
tually be measured in any experiment. The simplest, and possibly the most 
promising, type of experiments is deep-inelastic exclusive production of pho- 
tons, mesons, and perhaps even lepton pairs. Here we consider two experiments 
that have been studied extensively in the literature: deeply virtual Compton 
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scattering (DVCS) in which a real photon is produced, and diffractive meson 
production. There are practical advantages and disadvantages from both pro- 
cesses. Real photon production is, in a sense, cleaner but the cross section is 
reduced by an additional power of acm- The Bethe-Heitler contribution can be 
important but can actually be used to extract the DVCS amplitude. Meson 
production may be easier to detect, however, it has a twist suppression of l/Q^. 
In addition, the theoretical cross section depends on the unknown light-cone 
meson wave function. 



Deeply virtual Compton scattering was first proposed in Rcf. |186, 187 1 as a 
practical way to measure the off- forward distributions. Consider virtual photon 
scattering in which the momenta of the incoming (outgoing) photon and nucleon 
are q{q') and P{P'), respectively. The Compton amphtude is defined as 



P' 



p 



(130) 



where q — (q + q')/2. In the Bjorken limit, ~q^ and P • q — > cx) and their ratio 
remains finite, the scattering is dominated by the single quark process in which 
a quark absorbs the virtual photon, immediately radiates a real one, and falls 
back to form the recoiling nucleon. In the process, the initial and final photon 
helicities remain the same. The leading-order Compton amplitude is then 



1 



1 



9T i dx[ l_- + _^_ ^e^F,(a;,^,t,g2) 

' _i \x — ^ + le X + t, — le J ^ 



+ le^'^'^pa^np i dx ^ 

^1 \x — E^+ie x + E^—ie 



1 



1 



^e^^,(x,e,t,Q') 
(131) 



where n and p are the conjugate light-cone vectors defined according to the 
collinear direction of q and P, and is the metric tensor in transverse space. 
^ is related to the Bjorken variable Xb = —q^l{2P ■ q) by xb = 2^/(1 -I- 5). 

Much theoretical work has been devoted to DVCS in the last few y ears. The 
one-loop corrections to DVCS have been studied by Ji and Osborne [195|. An 



all-order proof of the DVCS factorization has been given in Rcf. |247, 248, 249 
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lOEL Suggestions have also been made to test the DVCS scatterin g me cha- 
nism |122|. Asy mme tries for polarized DVCS have been considered in |187| and 
reconsidered in 145 , |6^ . DVCS with double photon helicity flips have been 



investigated in Ref. [174, |67|. The estimates for cross sections have been made 
in Ref. |7|, |7|. 

Development on the experimental front is also promising. Recently, both 
ZEUS and HI collaborations have announced the first evidence for a DVCS 
signature |9^, and the HERMES collaboration has made a first measurement 
of the DVCS single-spin asymmetry |B0|. More experiments are planned for 



COMPASS, JLAB and other future facilities [ |121 |. 

Heavy quarkonium production was first studied by Ryskin as a way to mea- 
sure the unpolarized gluon distribution at small x [p56[. In the leading-order 
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diagram, the virtual photon fluctuates into a cc pair which subsequently scat- 
ters off the nucleon target through two-gluon exchange. In the process, the 
pair transfers a certain amount of its longitudinal momentum and reduces its 
invariant mass to that of a J/^*. The cross section is 

^(7* + P - J/* + P') = '-^4^^^aUQ')m,Q')f (132) 



dt 

where 



3acmQ^ 



(Q^ + M^)/4, M is the J/ip mass, and Fe+e- is the decay width 
into the lepton pair. The equation was derived in the kinematic limit s 3> 
^ Ajf^ ^ t and the Fermi motion of the quarks in the meson was neglected. 
Two other important approximations were used in the derivation. First, the 
contribution from the real part of the amplitude is neglected, which may be 
justifiable at small x. Second, the off-forward distributions are identified with 
the forward ones. 

The above result was extended to the case of light vector-meson production 
by Brodsky et al., who considered the effects of meson structure in perturbative 
QCD |89]. They found a similar cross section. 



da 
It 



{-f*N -> VN) 



47r^rvmval{Q)rjl (xg(x,Q^))^ 



where the dependence on the meson structure is in the parameter 



Vv = 7; 



dz 



4>^{z){ dzq^^z) 



(133) 



(134) 



zil-z) 

and (j)^ {z) is the leading- twist light-cone wave function. Evidently, the above 
formula reduces to Ryskin's in the heavy-quark limit {(j>^ (x) — 5{x — 1/2)). 

The amplitude for hard diffractive electroproduction can be calculated in 
terms of off- forward gluon distributions |247|. With the virtual photon and 
vector meson both polarized longitudinally (i.e. determined using a Rosenbluth 
separation, with the vector meson polarization measured via its decay products), 
one finds 



da 



dt 



(7*iV VN) 



2xf 



dx 



30^0111^ 
1 



1 



Fgix,^,t) 



(135) 



where again xb = 2^/(1 -f ^). The above formula is valid for any xb and t 
smaller than typical hadron mass scales. Hoodbhoy has also studied the effects 
of the off- forward distributions in the case of J/ip production |173|. He found 
that Ryskin's result needs to be modified in a similar way once the off-forward 
effects become important. 

More detailed theoretical studies of meson production have been done in 
Refs. |224, 16C, 149|. Longitudinal p° production data has been collected by 
the E665 and the HERMES collaborations and the comparison with model 
calculations is encouraging |274, 14S|. 
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8 Related Topics in Spin Structure 



In this section, we review two interesting topics related to the nucleon spin. 
First, we consider the Drell-Hearn-Gerasimov sum rule and its generalization to 
finite . Then we briefly review polarized A production from fragmentation of 
polarized partons where the A polarization can be measured through its weak 
non-leptonic decay. 



8.1 The Drell-Hearn Gerasimov Sum Rule and Its Gener- 
alizations 



The Drell-Hearn-Gerasimov (DHG) sum rule |126{ involves the spin-dependent 
photo-nucleon production cross section. Consider a polarized real photon of 
energy scattering from a longitudinally polarized nucleon and producing ar- 
bitrary hadronic final states. The total cross sections are denoted as era i(y\ 
where the subscripts 3/2 and 1/2 correspond to the helicity of the photon be- 
ing parallel or antiparallel to the spin of the nucleon. The sum rule relates 
the l/i/- weighted integral of the spin-dependent cross section from the inelastic 
threshold to infinity to the anomalous magnetic moment of the nucleon k. 



°° dv 27r2aemK^ 



-(.3(.)-a.(.)) = ^^. (136) 

For the proton and neutron, the sum rule is 204.5/i6 and 232.8//6, respectively. 

There has been much interest in recent years in testing the above sum rule by 
determining the integral on the left-hand side. Direct experimental data on the 
spin-dependent photoproduction cross section has become available recently p7| 
(see Fig. and more data at higher energy are coming soon However, 
many of the published "tests" in the literature rely on theoretical models for 
the photoproduction helicity amplitudes which are only partially constrained by 



unpolarized photoproduction data ||5^, 16£ | . Because of the Xjv weighting, the 



low energy amplitudes play a dominant role in the DHG integral |12S|. In fact, 
one can show that in the large limit, the integral is entirely dominated by 
the A resonance contribution |l07|. 



We will not discuss in detail how the phcnomcnological estimates of the DHG 



integral are done in the literature |201, 281 J94, 257 1. The interested reader can 



consult a recent review on the subject |127[ . The main conclusion from these 
calculations is that the isoscalar part of the sum rule (219/^6) is approximately 
satisfied, whereas a large discrepancy remains for the isovector part (— 14/x b). 
Typically, the proton integral is estimated to be in the range of 260/i6 to 290/Lt6. 



A more up-to-date analysis |12£] including the recent data from MAMI and 
the extrapolation of DIS data gives a result of 202 ± 10/i6 for the proton, but 
disagrees with the expected neutron sum by ^ 60/i6. 

What do we learn about nucleon spin physics by testing the sum rule? More- 
over, the DHG sum rule is the analogue of the Bjorken sum rule at Q2 = El 
(here, we discuss the Bjorken sum rule in the generalized sense that the first mo- 
ment of g\(x^ Q^) is related to nucleon axial charges in the asymptotic limit). If 
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Figure 23: Spin-dependent cross section for 7p 
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both sum rules are important to study, how do we extend these sum rules away 
from the kinematic limits {Q^ = and = oo)? Finally, how is the DHG sum 
rule evolved to the Bjorken sum rule and what can we learn from the evolu- 
tion? In recent years, there has been much discussion in the literature about the 
generalized DHG integrals and their dependence |71, 265, 130, 259 . A 

summar y of differ ent definitions of the gener alized DHG integrals can be found 
in Ref. |24C, 12£]. As pointed out in [197|, the key to addressing the above 
questions is the dispersion relation for the spin-dependent Compton amplitude 
5i(i.,Q2). 

The virtual-photon forward scattering tensor defines the spin-dependent am- 
plitude Si{i', uj), 

T"P{P,q) = i J e"'-^d^^{PS\TriOJ^iO)\PS) , 
where is the electromagnetic current. From general principles, such as causal- 
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ity and unitarity as well as assumptions about the laige-v behavior of Si{Q'^, v) 
one can write down a dispersion relation 



SMQ')-^r '^'^^f^ , (138) 



where G\{y' ^ Q^) is the spin-dependent structure function discussed in Sec. 1.3.1 
Whenever Si is known, in theory or experiment, the above relation yields a 
dispersive sum rule. For instance, the Bjorken and DHG sum rules are obtained 
from theoretical predictions for S'i(0,Q^) at = (X),0, respectively ^23[ . 

What do we learn by testing these dispersive sum rules? First, we learn about 
the assumptions required for the derivation of the relation; in particular, the 



high-energy behavior of the Compton amplitude |170|. Second, we learn about 
the scattering mechanisms in the virtual-Compton process. For the Bjorken 
sum rule, it is perturbative QCD and asymptotic freedom; for the DHG sum 
rule, it is nucleon-pole dominance and gauge symmetry p23| . Finally, if the 
sum rules are reliable, we have a new way to measure nucleon observables. In 
earlier sections, we discussed how to extract gA and AS (the fraction of the 
nucleon spin carried by quark spin) from polarized DIS data. Assuming the 
validity of the DHG sum rule, we obtain the magnetic moment of the nucleon 
from inclusive photoproduction. 

How do we extend these sum rules to other kinematic regions? According to 



Eq. |138| , the virtual Compton amplitude is the key. As discussed in Sec. |l.3| , at 
large but finite Q^, perturbative QCD introduces two types of corrections. The 
first are the radiative corrections: gluons are radiated and absorbed by active 
quarks, etc. The second are the higher twist corrections in which more than one 
parton from the target participates in the scattering. With these corrections, 
we can extend Bjorken's result for the Compton amplitude from = oo to 



finite ]194 193 1. Since the scale that controls the twist expansion is on 
the order of 0.1 — 0.2 GeV^, the perturbative QCD prediction for S'i(0, Q^) is 
valid down to ~ 0.5 GeV^. Combined with Eq. 138, it yields a generalized 



Bjorken sum rule. It is the generalized Bjorken sum rule that is commonly 
tested experimentally. 

At small but finite Q^, chiral perturbation theory provides a sound the- 
oretical method to calculate corrections to the low-energy theorem ||7l|, 197| 



Recently, a fourth-order chiral perturbation theory calculation for the inelastic 



part of Si{Q,Q^) yielded |191| 



5i(0, Q') = -k" + .^./^/{a (1 + 3Ky + 2(1 + iKs)T-')Q' + (139) 

12(47r/jr)''m^ 

The result shows a rapid Q^-dependence near ~ 0, which is qualitatively, 
though not quantitatively, consistent with a recent phenomenological analy- 



sis |129|. For a quantitative test, one needs polarized electron scattering data 

soon available from JLab ]95[ |. 

How does the DHG sum rule at = evolve to the Bjorken sum rule at 
= oo? The physically most interesting quantity which connects both sum 
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rules is 



where Fi{Q'^) and F2{Q^) are the elastic nucleon form factors. It is the elastic 



contribution which dominates at low [197|. r((5^) starts at 1(0) + k/2 from 



the proton (neutron) at = and rapidly decreases to about 0.2 at — 0.7 
GeV^ and remains essentially flat as oo. The interpretation for the 



variation is as follows 1 197 1 . The forward Compton amplitude is an amplitude for 
the photon to scatter from a nucleon target and remain in the forward direction. 
This is very much like a diffraction process and T{Q^) is the "brightness" of the 
diffraction center. For low photons, scattering from the different parts of 
the proton is coherent, and the scattered photons produce a large diffraction 
peak at the center. As becomes larger, the photon sees some large scale 
fluctuations in the nucleon; the scattering becomes less coherent. The large 
scale fluctuations can largely be understood in terms of the dissociation of the 
nucleon into virtual hadrons. When > 0.5 GeV^, the photons see parton 
fluctuations at the scale of l/Q. As ^ oo, the photons see individual quarks 
inside the nucleon and the scattering is completely incoherent. The diffraction 
peak is just the sum of diffractions from individual quarks. In short, the 
variation of the sum rules reflects the change of the diffraction intensity of the 
virtual photon as its mass is varied. 

A clear theoretical understanding of the virtual photon diffraction at ~ 
0.1 — 0.5 GeV^ is not yet available, but there are two distinct possibilities. 
First, there is a gap in which neither parton nor hadron language describes 
the scattering well. In this case, an interesting theoretical question is how 
the transition from low to high happens. Second, some extensions of the 
twist expansion and chiral perturbation theory may overlap in the intermediate 
region. If so, we have parton-hadron duality at a new level. In any case, a 



lattice calculation of 6*1(0, Q ) may shed important light on this [|19i 



8.2 Spin-Dependent A Fragmentation 

In the constituent quark model, the spin structure of the A baryon is simple: 
the ud quark pair couples to give zero angular momentum and isospin, and the 
spin of the A is entirely carried by the spin of the remaining s quark. From 
our present knowledge of the spin structure of the nucleon, we expect that this 
naive picture will fail to explain the actual spin structure of the A. In fact, 
if SU (3) flavor symmetry is valid, we can deduce from the beta decay data 
and polarized deep-inelastic scattering on the nucleon that (in leading order) 

AwA = AdA 0.23 and Asa ~ 0.58 [||. 

Unfortunately the spin structure of the A cannot be measured because of the 
lack of a stable target. However, the spin-dependent fragmentation of partons 
to the A baryon can be studied experimentally because the A polarization can 
be measured through the self-analyzing decay A . The fragmentation 

functions are difficult to calculate in QCD, even in principle. We have little 
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experience in modeling the fragmentation functions compared with the internal 
structure of the nucleon. Nevetherless, one hopes that the spin physics in the 
fragmentation process corroborates what we learn about the spin structure. 
Moreover, if a A or A is exclusively produced from the fragmentation of a strange 
or antistrange quark, respectively, the measurement of the A polarization is a 
way to access the strange quark polarization in the nucleon. 

A relatively simple process from which the spin-dependent fragmentations 
to A can be studied is e~^e~ annihilation with one of the beams (say, electron) 
polarized. Considering only the intermediate photon state, the asymmetry in 
polarized A production is 



-^cos0j24{^D,iz) + AD,{z)) , (141) 



d^Aaie-e+ AX) _ al„ 

dildz 2s ^ 

where A.Dq{z) = D'^{z)~D'^ (z) is a spin-dependent fragmentation function and 
D^{z) are the fragmentations of the quarks with helicities ±1/2 to a A of helicity 
+1/2. At the peak, the parity violating coupling induces polarizations in 
the quark-antiquark pairs produced. Hence even without beam polarization, 
the A particles produced through fragmentations are polarized jo^. Recently, 
several collaborations at LEP have extracted the A polarization from quark 
fragmentation at the Zq peak [ p7| . A number of models for spin-dependent 



quark fragmentation functions have been proposed to explain the results [208 



116, qO| 22£], and data are consistent with very different scenarios about the 
flavor structure of fragmentation. 

The polarized fragmentation functions can also be measured in deep-inelastic 
scattering, in which the polarized beam produces a polarized quark from an 



unpolarized target, which then fragments |178|. Within the QPM, the measured 



A polarization from a lepton beam with polarization Pi is, 

where D{y) is the depolarization factor. A process-independent A polarization 
can be defined from Pa = Pcxp/ {PbD{y)). 

The A polarization from DIS scattering was first measured by the E665 
Collaboration with a 470 GeV/c^ polarized muon beam (F^ = -0.7 ± 0.1) ||l||. 
The data sample was taken at 10""' < xb < 10~^ with (xb) — 5 • 10~^, 0.25 < 
Q2 < 2.5 GeV^ with (Q^) = 1.3 GeV^, and {1^) = 150 GeV. The A polarization 
was found to be -1.2 ±0.5 at < < 0.3 and -0.32 ±0.7 at 0.3 < xp < 1.0. 
The A polarization was 0.26±0.6 and 1.1±0.8 for the two bins, respectively. The 



comparisons with different fragmentation models can be found in Refs. [£6, 230 1. 

Recently, HERMES has also reported a measurement of the A polarization 
from polarized deep-inelastic positron scattering from an unpolarized proton 
target. The result is Pa = 0.11 ± 0.17 ± 0.03 at an average z = 0.45 |l|- The 
result seems to be consistent with the assumption of the naive quark model that 



the A polarization is entirely carried by the valence s quark \llt\ 
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In Ref. 1 117 1, predictions for A production from p — p collisions at RHIC and 
HERA-iV with a single beam polarization was studied. Spin asymmetry mea- 
surements as a function of the rapidity provide a way to discriminate various 
models of the spin-dependent fragmentation. The main theoretical uncertain- 
ties, such as the NLO corrections and the unknown polarized parton distribu- 
tions, have no major impact on the asymmetry. In Ref. [8^ , it is argued that the 
hyper-fine interaction responsible for the N-A mass splitting induces a sizable 
fragmentation of polarized up and down quarks into a A, which leads to large 
positive A polarizations at large rapidity. 

9 Conclusions 

Since the EMC publication of the measurement on the fraction of the nucleon 
spin carried by quarks, understanding the spin structure of the nucleon has 
become an important subfield in hadron physics. In this review, we have tried 
to highlight some of the important developments over the last ten years and 
discuss some of the future prospects in this field. 
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